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The Gaia mission is designed as a Galaxy explorer, and will measure simultaneously, 
in a survey mode, the five or six phase space parameters of all stars brighter than 
20th magnitude, as well as providing a description of their astrophysical characteris- 
tics. These measurements are obtained by combining an astrometric instrument with 
micro-arcsecond capabilities, a photometric system giving the magnitudes and colours 
in 15 bands and a medium resolution spectrograph named the Radial Velocity Spec- 
trometer (RVS). The latter instrument will produce spectra in the 848 to 874 nm 
wavelength range, with a resolving power R = 11500, from which radial velocities, 
rotational velo cities, atmospher ic parameters and abundances can be derived. A com- 
panion paper ijKatz et alJl2004() presents the characteristics of the RVS and its per- 
formance. This paper details the outstanding scientific impact of this important part 
of the Gaia satellite on some key open questions in present day astrophysics. The 
unbiased and simultaneous acquisition of multi-epoch radial velocities and individual 
abundances of key elements in parallel with the astrometric parameters is essential 
for the determination of the dynamical state and formation history of our Galaxy. 
Moreover, for stars brighter than V ~ 15, the resolving power of the RVS will give 
information about most of the effects which influence the position of a star in the 
Hertzsprung-Russcll diagram, placing unprecedented constraints on the age, internal 
structure and evolution of stars of all types. Finally, the RVS multi-epoch observations 
are ideally suited to the identification, classification and characterisation of the many 
types of double, multiple and variable stars. 

Key words: binaries: spectroscopy - stars: general - Galaxy: evolution - Galaxy: 
formation - Galaxy: kinematics and dynamics - Galaxy: structure 
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1 INTRODUCTION 

The Gaia mission was approved as one of the next corner- 
stones of the European Space Agency (ESA) scientific pro- 
gramme in September 2000, and was confirmed as one of the 
flagship missions of ESA C osmic Vision in November 20 03. 
The mission was conceived dLindeeren fc Perrvmanll996l) in 
the light of the success of Hipparcos, the first space astrome- 
try mission, also operated by ESA. This first mission proved 
the very high potential of such an instrument by providing, 
for the first time, a large data set of high-accuracy, homoge- 
neous absolute trigonometric parallaxes and proper motions. 
Milli-arcsecond level accuracy was obtained for 118 000 pre- 
selected stars. More than 3000 papers using this unique set 
of data have been published to date, in domains ranging 
from reference frames to the cosmic distance scale, through 
stellar physics and Galactic astronomy. 

Three aspects of the Gaia mission represent major im- 
provements with respect to Hipparcos: an unprecedented as- 
trometric accuracy of a few micro-arcseconds, the system- 
atic, and hence unbiased, observation of some one billion 
objects down to a G magnitude 1 of 20, and the simulta- 
neous acquisition of astrometric parameters and of multi- 
epoch multi-colour photometry and spectroscopic observa- 
tions. The overal l science case for Gaia is impressive and 
multi -faceted (see lPerrvman et al l200lHBienav me fc TurorJ 
120021 . and references therein). It will clarify the origin and 
formation history of our Galaxy through the provision of 
a quantitative stellar census, a kinematical and a dynami- 
cal description of each of the Galactic stellar populations, 
the identification of accretion debris, a determination of the 
Galactic rotation curve and the distribution of dark mat- 
ter. It will also detect and characterize all types of stars, 
even the rarest types or the more rapid evolutionary phases, 
placing unprecedented constraints on the internal physics 
of stars, stellar evolution and stellar ages. The Gaia results 
will also firmly establish the distance scales in our Galaxy 
and far beyond and will detect tens of thousands of double 
and multiple systems, variable stars, minor bodies in the 
Solar System, quasars and supernovae. Thousands of extra- 
solar planets will be systematically identified and their pre- 
cise orbits and masses determined. Finally, Gaia will per- 
form a strong test of General Relativity by measuring the 
Parametrized Post-Newtonian (PPN) parameters 7 and j3 
with unrivalled precision. 

The Gaia spectroscopic instrument, the Radial Veloc- 
ity Spectromete r (RVS), will co ntribute significantly to this 
science case fsee lMunarill2003bl . and references therein), by 
providing spectra for 100 to 250 million stars up to magni- 
tude V = 17. Its key contributions will be the following: 

• the simultaneous acquisition of the third component of 
the space velocity, the radial velocity, essential for kinematic 
and dynamical studies in the Galaxy. These data will give 
direct insight into the past accretion events which led to the 
formation of the Galaxy as it is now observed, and will place 
strong observational constraints on formation scenarios for 
the various components of the Galaxy; 



1 The characteri stics of the G filter are discussed in detail in 
iKatz et al ] l2004t hereinafter Paper I). For example, a G5V star 
has G — V = -0.36. 



• multi-epoch measurements of radial velocities, essential 
for the detection, classification and complete characterisa- 
tion of binary systems, and key observables for the determi- 
nation of the fundamental parameters of pulsating variable 
stars; 

• the determination of the 'perspective acceleration', 
mandatory for obtaining the expected astrometric accuracy 
for high velocity stars; 

• the astrophysical characterisation of the brightest sub- 
set of Gaia stars, complementary to the information pro- 
vided by Gaia photometric observations. In particular, the 
overall metal abundance will be derived up to V ~ 14-15 
and individual element abundances will be obtained up to 
V ~ 12-13. For fainter stars, photometry and astrometry 
will be used to estimate abundance parameters - RVS spec- 
tra will calibrate these estimates at the bright end of the 
luminosity function. These data will help to resolve the age- 
metallicity degeneracy which currently hinders interpreta- 
tion of the colour-magnitude diagrams that are a fundamen- 
tal input to our understanding of the chemical evolution of 
the Galaxy; 

• the observation of the 862 nm Diffuse Interstellar 
Band which will contribute to the construction of a three- 
dimensional map of Galactic interstellar reddening. 

The inclusion of the RVS on board the Gaia satellite 
has critical advantages over any ground-based survey be- 
cause it will obtain spectra simultaneously with the astro- 
metric parameters and photometric indices, and for system- 
atically the same stars. Even though much effort was put 
into the planning of spectroscopic ob servations of a large 
proportion of the H ipparcos target list iGerbaldi et al .119891 : 
iMavor et all I1989T) with particular emphasis on the kine- 
matically unbiased selection of the stars to be observed, the 
large number of required observations has meant that the 
resulting measurement programme is still incomplete and 
only parts of it have be en published dGrenier et alJll999bO 
iNordstrom et alj|2004 - 

The absence of an un-biased radial 
velocity data set to complement the Hipparcos astrometric 
observations constituted a severe constraint on the value of 
th e data set for the st udy of Galactic kinematics. In fact, 
as iBinnev et al] <ll997t) emphasised, the non-uniformity of 
published velocities for stars in the Hipparcos sample meant 
that (in 1997) the set of stars with full three-dimensional 
velocity information was kinematically biased because high 
proper motion objects were more likely to have measured 
radial velocities. This had a significant impact on the poten- 
tial scientific value of this data set. In the case of Gaia, the 
measurement principles ensure the total homogeneity of the 
observations, with the sole exception of the approximately 
10 percent of the densest regions of the Galaxy (the cores 
of globular clusters and the inner Galactic bulge) . After the 
detection of particularly interesting objects, it will clearly 
be fruitful to observe them with ground-based instruments 
providing much higher spectral resolution. 

This paper summarises the expected impact of the 
RVS on our knowledge of Galactic structure and evolution 
(§Hj; binary star statistics and characterisation (§|3Jl; stel- 
lar physics and evolution, as well as variable stars (§|1J. In 
§ l2.1l we briefly summarise the performance of the RVS in es- 
timating velocities as a function of position on the sky as this 
determines the contribution which the RVS will make to our 
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Table 1. Precision (la error) of the RVS instrument in the de- 
termination of stellar parameters. The table gives the expected 
accuracies for effective temperature T e g, surface gravity log g, iron 
abundance [Fc/H], and rotational velocity vsini. For all param- 
eters other than rotational velocity, the quoted precisions will be 
obtained for a G5V star to a limiting magnitude of V~ 14.0. 
Limiting magnitudes for other stellar types, as well as a complete 
discussion of the algorithms used to calculate these accuracy esti- 
mates, will be presented elsewhere (Recio-Blanco et al., in prepa- 
ration) . Paper I also contains some additional discussion of these 
performances. 



Expected stellar parameter accuracies 

a(T eS ) = 75 ± 55 K 
cr(logg) = 0.16 ± 0.12 dex 
o-QFe/H]) = 0.11 ± 0.08 dex 

a(vs'mi) = 5kms _1 to V=15 for late type stars 

cr(tJsini) = 10 — 20kms~ 1 to V=10-ll for B5 main-sequence stars 



Table 2. Precision (la error) of the RVS instrument in the de- 
termination of radial velocities. The quantity a\ is the la error in 
radial velocity at the limiting magnitude V\ after a single transit 
of the RVS instrument; cr2 is the la error in radial velocity at 
the limiting magnitude V2 after 102 transits (= average number 
of transits per star during the mission); V3 is the limiting mag- 
nitude for velocity errors of lkms~ 1 after 102 transits. Results 
for four stellar types are shown: an F2 II, a G5V and a Kl III 
star each of solar metallicity and a low-metallicity Kl III star. In 
each case, the limiting magnitude V2 is taken to be the faintest 
magnitude (to the nearest 0.5 magnitudes) at which the radial 
velocity error is < 20kms" 1 . For more details see Paper I. 



Expected radial velocity accuracies 



Type 


[Fc/H] 


v sin i 


ai 


Vi 


a 2 


v 2 


v 3 


F2 II 


0.0 


20 


14.4 


13.5 


20.8 


16.5 


13.0 


G5V 


0.0 


5 


19.7 


14.5 


16.0 


17.0 


14.0 


Kl III 


0.0 


5 


11.5 


15.0 


10.6 


17.5 


14.0 


Kl III 


-1.5 


5 


11.7 


14.5 


10.3 


17.0 


14.0 



understanding of the various components of the Galaxy. Ta- 
bles0and|5|summarise the expected precisions and limiting 
magnitudes for various stellar parameters including radial 
velocity. For more details of the performance of the RVS in- 
strument, the reader is directed to Paper I. In the present 
paper, we emphasise the most important potential scientific 
contributions of the RVS. However, we note that there will 
be other topics, in addition to those discussed here, for which 
the RVS data set will be valuable. 



2 GALACTIC STRUCTURE AND EVOLUTION 

One of the key science drivers for the Gaia mission is to de- 
termine the current dynamical state and formation history 
of the Milky Way. Radial velocities provide the 6th phase 
space coordinate, vital for building a complete picture of 
the stellar dynamics in the Galaxy. In this regard, an unbi- 
ased sample of radial velocities is an essential complement 



to the Gaia proper motion and distance data set. The pres- 
ence of the RVS on board the Gaia satellite will ensure that 
multi-epoch radial velocities will be acquired in an un-biased 
manner for all stars brighter than V ~ 17. 

The most popular scenario for the formation of the 
Milky Way involves the merger of several similar sized sys- 
tems (whose exact nu mber may vary between 5 and 20) up 
to a redshift 2 ~ 2 — 3 llKauffmann et al ] ll993tlAbadi et al] 
120031) . These systems contribute both dark matter, stars and 
gas. The accretion of smaller systems occurs throughout the 
formation history of the Galaxy, but it must have been the 
dominant form of mass growth since z ~ 1, when the thin 
disc probably began to form. It is unclear whether the bulge 
of the Galaxy was formed from the mergers of equal size ob- 
jects at very high z, or from the instability of a pre-existing 
disc. At high z, when mergers were more frequent, it seems 
quite likely that discs could easily become unstable, lead- 
ing to the formation of a bar/bulge. The mass accreted at 
later times builds up t he external halo and the thin disc 
jGilmore fc Wvseil200ll) . In recent years, considerable evi- 
dence of past accretion events has been obtained (see S I2.2H . 
However, numerical simulations suggest that this scenario 
for the formation of the Galaxy is not without its prob- 
lems. For example, there is an apparent discrepancy between 
the hundreds of smaller sub-haloes surrounding a simulated 
Milky Way-type galaxy in a standard ACDM universe and 
the 11 dwarf satelli te galaxies which a r e observed to orbit 
the Milky Way (e.g. lMoore et alJll999l) . lAbadi et al] <2003l) 
also raise the concern that simulated ACDM disc galaxies 
typically have significant bulge components which may be 
difficult to reconcile with the observed numbers of nearby, 
pure-disc galaxies. It is therefore important to place observa- 
tional constraints on the possible formation histories of the 
various Galactic components. The most direct way to do this 
is to look for correlations between the ages, metallicities and 
space motions of their associated stellar populations. 

In this section, we highlight some of the key Galactic 
stellar populations which will be amenable to study with the 
RVS and the scientific issues which can be addressed using 
the RVS data set. We note that this is not intended to be an 
exhaustive review of Galactic structure or of the contribu- 
tion of the Gaia mission as a whole to our understanding of 
the Galaxy, but rather a discussion of the key contributions 
in this area to be expected from the RVS instrument. 

2.1 Expected RVS velocity performance 

Fig. Qshows the stellar number density for all stars brighter 
than F= 17.5 as a function of Galactic coordinates I 
and b taken from th e Guide Star Catalogue II (GSCII: 
iBucciarelli et alj|200ll) . This spectral band is close to that 
of the RVS and in what follows we assume that these num- 
ber densities are representative of those which the RVS will 
encounter. As is discussed in Paper I, crowding at low Galac- 
tic latitudes severely limits the ability of the RVS to provide 
robust vel ocity measurem ents. Simulations of the effects of 
crowding ]Zwitterl2003bl) have demonstrated that degrada- 
tion of radial velocity determinations sets in when the num- 
ber density exceeds about 2 x 10 4 stars per square degree. 
In the Figure, red and black correspond to stellar densities 
above 2 x 10 4 stars per square degree and 4 x 10 4 stars per 
square degree, respectively. Thus, it is likely that the region 
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Figure 1. Stellar density as a function of Galactic coordinates (l,b) for stars brighter than F= 17.5 in the GSCII catalogue. Note that 
the typical V— F colour of the disc stellar populations is about 0.6. Different colours correspond to regions in which the stellar density 
is (i) less than 5 X 10 3 deg -2 (white) (ii) 5 X 10 3 — 2 X 10 4 deg~ 2 (green) (hi) 2 - 4 X 10 4 deg~ 2 (red) and (iv) above 4 X 10 4 deg -2 
(black). Current simulations show that the blue and green regions will be fully observable by the RVS. However, the black regions will 
be inaccessible as well as some of the red regions. 





Figure 2. Expected limiting magnitude Fi; m for RVS radial velocities as a function of Galactic coordinates (I, b). The colours correspond 
to limiting magnitudes of F lim = 17.5 (blue), F lim = 17.25 (green), F lim = 16.75 (red), F lim = 16.25 (light grey), F Um = 15.25 (dark 
grey) and Fi; m = 14.25 (black). Note that this plot does not take into account telemetry issues which will affect the densest regions of 
the Galaxy. See text for discussion. 



of the sky within 10 degrees of the Galactic centre will not 
be observable with the RVS. Techniques for the extraction 
of radial velocities from the data in crowded fields are still 
being investigated, with the prospects of some (limited) im- 
provements over what is currently achievable. We also note 
that the stellar densities in the Figure are av eraged over tiles 
of on e square degree. It is well known (e.g. Popowsk i et al.l 
I200&1) that extinction levels towards the Galactic centre vary 
on small scales, and so regions of high stellar density may be 
adjacent to regions of significantly lower density. However, 
since the lowest extinction regions (e.g. Baade's window) 
have stellar densities exceeding 10 6 stars per square degree, 
Fig. Q probably gives a realistic reflection of the regions in 
which the RVS will be able to observe. 

For stars of a given apparent magnitude Fn m , simula- 
tions have shown that, for late type stars, severe degradation 
of the radial velocity accuracy occurs if the total stellar den- 
sity of stars to a limit roughly 0.75 magnitude fainter than 
Fii m exceeds 2 x 10 4 stars per square degree. Thus, for a given 



region of the sky, it is possible to determine an effective mag- 
nitude limit above which we can expect the RVS to provide 
accurate radial velocities. Fig.|2]shows the effective limiting 
magnitude for radial velocity measurements as a function of 
position on the sky. This shows that over most of the sky the 
RVS will be able to obtain velocities to a magnitude limit 
of between Fn m = 16 and Fn m = 17, although in the direc- 
tion of the Galactic centre this limit is significantly brighter. 
Given that the disc populations have a mean V — F colour 
of 0.6, the limiting V magnitude will range between 16.85 
to 17.5 over most of the sky. In fact, telemetric constraints 
will severely restrict the amount of data from the Galactic 
centre region which can be obtained and so the magnitude 
limit for these regions in Fig. [5]is probably optimistic. In the 
following, we assume that the RVS performances described 
in Paper I will be achieved over approximately 90 per cent 
of the sky, as suggested by Figs. [5] and 
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2.2 The Stellar Halo and Debris Streams 

If the Milky Way was built via accretion as in t he hier- 
archical formation model (e.g. IWhite fc Reeslll978fl . accre- 
tion events should have left fossil signatures in its present 
day c omponents, which should be clearly dete ctable with 
Gaia iHelmi fc Whitdll999l lAbadi et al.fl2003l) . The most 
natural place to look for such substructures is the Galac- 
tic stellar halo, since a spheroidal component is formed 
by the trails of stars left by disrupted sat ellite galax- 
ies l|.Tohnston et alJll99.4 iMorrison et alJl2000t ). Moreover, 
it is here where the most metal-poor stars are found, which 
cou ld imply that this component was in place very early 
on ijFreeman fc Bland-HawthorrjEio3) . 

Recent observations have shown that indeed consider- 
able structure is still present in the Milky Way's halo, in- 
dicating that accretion events have had s ome role in its 
formation history (e.g. W94 IMaiewski et alJ 



19961: iHelmi et al.lllQ9fll: llvezic et alJl2000l) . The existence 
of a substantial counter-rotating halo component is sug- 
gested by RR- Lyrae and blue horizontal b ranch (BHB) 
star kinematics l|Maiewskilll992l : ICarnevlll999l) . In the Solar 
neighbourhood, it is unclear whether the local halo veloc- 
ity d istribution is well approximated by a smooth Gaus- 
sian ijMartin fc Morrisonl Il99l IChiba fc Beersll200(tl which 
might point to the presence of substructure. A large frac- 
tion of the substructures discovered so far are located in the 
outer regions of our Galaxy, where the debris can remain 
spatially coherent for many Gyr. This implies that it can 
be recovered from low dimensionality surveys, such as 2D 
(using two sky coordinates), 3D (using also distance) or 4D 
(including radial velocities) maps. Many surveys have been 
able to discover debris using this approach, which requires 
suitable tracers such as giant stars (BHB, RR Lyrae, red gi- 
ant branch (RGB)) to probe the outer Galaxy. These surveys 
inclu de the Sloan Digital Sky Survey (SP SS: [ Yanny et al" 
20^ ), the Spaghetti Project Survey (SPS: IMorrison et a 
2000li . the Qu asar Equatorial S urvey Team (QUEST) RR 



Lyrae survey dVivas et al.1 1200 if), the Tw o-Micron All-Sky 
Survey f2MASS: lR.ocha-Pinto et al.l2003h and a few smaller 
ones. 

One of the most spectacular findings has been the dis- 
covery of debris from the Sagittarius d warf galaxy, which 
spans a complete great circle on the sky iRocha-Pinto et"all 
120031) ■ We are presently witnessing the merging of a 
satellite galaxy, which is just now contributing stars 
to build up the outer stellar halo of our Galaxy. In 
fact, most of the halo substructure discovered so far 
can be attributed to the Sg r dwarf l|lbata ct al. 2001b: 
Martmez-Dekado et all 1200 ll: iDohm-Palmer et alJ 12001 



Newberg et al.l2002l:lMarti'nez-Delgado et all2004T) . It is un- 
clear whether other accretion events in the outer halo still 
await discovery, although this should probably be the case 
if the Galaxy was built hierarchically. The present situation 
may be understood from the fact that the streams discovered 
so far are dynamically the youngest, and hence the easiest to 
detect. It is quite likely that other debris has escaped detec- 
tion because of its low surface brightness. Observationally, 
we need to keep in mind that no surveys have yet reached 
both the required depth and sky coverage. An ideal survey 
would need to have a magnitude limit not brighter than V 
= 19 or 20, excellent photometry to allow identification and 



reliable distance estimation, preferably spectroscopic follow 
up to measure radial velocities to an accuracy of 10 kms , 
and its sky coverage should allow the tracing of tidal trails 
over large portions of the sky, which would require at least 
a few thousand square degrees. 

Most of the stars in the Galactic halo are however, in- 
side the Solar circle. Thus, it is the inner halo that contains 
most of the information related to the merging history of 
our Galaxy. However, in this region of the Galaxy, the de- 
bris from a disrupted satellite galaxy phase-mixes rapidly 
(a simple consequence of the shorter orbital time scales and 
the more flattened nature of the Galactic potential), eras- 
ing s patial correlations almost completely iHelmi 



White! 



119991) . Nevertheless, the debris still retains its coherence in 
phase-space which is reflected in its kinematics. For example, 
if the whole stellar halo was built from merged satellites, we 
would expect between 300 and 500 (mainly) cold streams 
in the Solar neighbourh ood, whose origin could be traced 
back to those satellites iHelmi fc Whitdll999l : IHelmi et alJ 
l2003al) . 

To test this prediction, large samples of halo stars with 
accurate 3D velocities are required. For example, to resolve 
the stellar halo velocity ellipsoid near the Sun into the in- 
dividual streams at a 3-<r level, implies that the required 
accuracy in any velocity component e v should satisfy 



3™ N * 

O iv s tr 



(1) 



Here owo is the typical velocity dispersion of halo stars near 
the Sun and is approximately 100 kms -1 , A^treams is the 
number of streams predicted in a small volume around the 
Sun, which is of order 500, and n is the dimensionality of the 
kinematic survey. Thus, if radial velocities alone are avail- 
able, then the required accuracy is e v < 0.07 km s - . If only 
the tangential velocities are available, then e v < 1.5 kms - , 
while in the case where all 3 components are obtainable, 
e v < 4.2 kms -1 , which is within the reach of the RVS (for a 
K1III star, this precision will be achievable for stars brighter 
than V = 16.0 — 16.5 depending on metallicity) . Clearly, 
knowledge of the radial velocities is critical to discover the 
streams and measure their properties. It is worth empha- 
sising that kinematically cold streams are expected even in 
the full hierarchical regime, where mergers of systems of 
more equa l mass, rather than simple satellite accretion, are 
dominant (IHelmi et aj2003rj) 2 . The dumpiness in the kine- 
matics of halo stars should thus be a distinct feature of the 
hierarchical formation of our Galaxy. 

Since an individual disrupted satellite may give rise to 
many streams in a small volume around the Sun (for exam- 
ple, an object of the size of the Small Magellanic Cloud could 
produce about 30 str eams after orbiting the inner Galaxy for 
a Hubble time, fsee lHelmi fc Whitdll999T) l. more sophisti- 
cate d methods are needed to id entify the debris observation- 
ally. IHelmi fc de Zeeuwl J^OOPi) proposed a method based on 
the expected dumpiness in the space of the integrals of mo- 
tion (energy, angular momentum) for stars having a common 



2 This is a consequence of the conservation of phase-space den- 
sity: as tidal debris spreads out in coordinate (physical) space, the 
local velocity dispersion decreases, so cold streams are expected 
if the mergers took place a relatively long time ago. 
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origin. This method has been successfully applied to the So- 
lar neighbourhood to identify an ancient minor merger that 
contr ibuted 10 per cent of the nearby halo stars ijHelmi et alJ 
Il999h. 

iHelmi fc de Zeeuwl feOOd) simulated the entire stellar 
halo of the Galaxy starting from disrupted satellite galax- 
ies 3 , and "observed" it with the Gaia mission by convolving 
the positions and velocities of the particles in the simula- 
tions with the expected errors. They showed that the initial 
clumping in the space defined by energy E, total angular 
momentum L and its ^-component L z , is maintained to a 
great extent even after 12 Gyr of evolution and the error 
convolution. Fig. |2] presents the results of observations of 
the same simulations but incorporating the most recent er- 
ror estimates for the observed parallaxes and the expected 
RVS radial velocity errors. In a ddition, a magnitude limit 
of V = 17.5 has been assumed iHelmi fc de Zeeuwl (1200(f) 
assumed a limit of V — 15). 

Fig. [3] illustrates the importance of radial velocities in 
the identification of the streams in the Solar neighbourhood. 
The top left-most panel in the Figure shows the sky distri- 
bution of particles from a simulation of a disrupted satel- 
lite orbiting for 8 Gyr in the outskirts of the Galaxy with 
an apogalactic distance of about 60 kpc and a perigalac- 
ticon of about 20 kpc. The particles correspond to red gi- 
ant stars (each assumed to have an absolute magnitude of 
Mv = +1) located within 40 kpc of the Sun (after error 
convolution). The region highlighted in blue contains only 
those stars which are within 10 kpc of the Sun. The panel 
beneath shows their motions (v a ,vs) as would be derived 
only on the basis of their proper motions and parallaxes. It 
is clear that the astrometric information alone is useful for 
finding debris in the outer halo. 

The middle panels in Fig. [3] illustrate the case of a 
stream in the inner halo. The top panel shows the sky distri- 
bution of particles from a simulation of a satellite orbiting in 
the inner Galaxy with an apogalacticon of about 12 kpc and 
a perigalacticon of about 4 kpc. The bottom panel shows the 
(v a ,vs) distribution of the stars located in a sphere of 1.5 
kpc radius around the Sun and for which the velocity errors 
are less than 25kms~ 1 and relative distance errors less than 
30 per cent. Clearly, astrometric information alone is not 
sufficient to find debris streams in the inner halo. 

Given that most of the stellar halo mass is actually lo- 
cated in the inner halo, it is vital that we are able to identify 
stellar streams in this volume also. The top right panel of 
Fig- ED shows the distribution of energy E and angular mo- 
mentum about the z-axis L z for all the stars from the inner 
halo s atellite after error convolution. As IHelmi fc de Zeeuwl 
teOOOf) found, the stars retain their clustering in the space 
defined by the integrals of motion. For each star in the RVS 
sample, both E and L z can be computed. The bottom right 
panel shows the (L z , E) distribution for stars originating in 
33 different satellites that could be populating the present 



3 Their simulations of satellite disruption adopted a fixed Galac- 
tic potential, and hence are quite idealistic and may not be very 
represent ative of the c o nditio ns found in a hierarchical universe. 
However, iKnebe et alJ i2005T) used cold dark matter simulations 
with a live Galactic potential and found that substructures do 
remain quite coherent in energy-angular momentum space. 



day stellar halo. For each star E and L z have been com- 
puted after proper error convolution on all observed quanti- 
ties. A standard friends-of-friends algorithm applied to the 
(L z , L, E) space (where L is the magnitude of the total an- 
gular momentum vector) is able to recover two thirds of all 
the satellites. The use of 6-D information is critical for the 
recovery of these accretion events: this would not be possi- 
ble without the radial velocity information provided by the 
RVS. 

More realistic, cosmologically motivated, models of the 
Galaxy and in particular of its stellar halo, as well as a 
deeper understanding of the phase-space structure of the 
debris from mergers, need to be produced. These would en- 
able us to develop more refined methods to recover the de- 
bris observationally and, once the Gaia data are available, 
to unravel the formation history of our Galaxy. 

2.3 Dark Matter in the Milky Way 

Our knowledge of the total mass and extent of the dark 
halo of the Milky Way is sadly lacking given its impor- 
tance for compa r isons with theories of galax y formation (e.g. 
iKochanekl Il99r3 : IWilkinson fc Evanl Il999h . The RVS will 
significantly enhance our understanding of the dark matter 
distribution within 50 kpc of the Galactic centre. 

The volume and column mass density of the Galactic 
disc place limits on the nature of the dark matter which 
makes up the halo. Locally, we can place robust limits on the 
mass in stars and gas based on star counts and our knowl- 
edge of the stellar mass function. The total mass implied by 
models of the stellar kinematics can be compared with the 
observed mass f ee. iKuiiken fc GilmorJ ll98St I Creze et al.l 
ll998tlChen et al.ll2003ti . Current estimates suggest that the 
two approaches yield similar values, which argues against 
the presence of a significant amount of dark matter in the 
Galactic plane and therefore the existence of dissipational 
forms of dark matter. The scale-height of any unseen com- 
ponent can be inferred from observations of populations at 
a range of radii and heights above the Galactic plane which 
wil l show the variation of vo lume density within the disc. 
As lKuiiken fc Gilmord (I1989T) discuss, a key uncertainty in 
current estimates of the local mass density lies in the uncer- 
tain variation in the "tilt" of the stellar velocity ellipsoid as 
a function of height or, in other words, whether the velocity 
ellipsoid is aligned with the axes of a spherical (maximal 
tilt) or cylindrical (no tilt) coordinate frame. The accuracy 
of the RVS velocities will be sufficient to map out the veloc- 
ity structure over a large volume surrounding the Sun. For 
example, the radial velocities of metal poor K1III stars will 
be measurable with accuracies of < 5kms _1 to a limiting 
magnitude of V ~ 16, corresponding to a distance of 10 kpc 
(assuming My = 1). Radial velocities are necessary to de- 
termine the detailed shape of the velocity distribution in a 
model-independent way. 

Another important quantity related to the mass of the 
Milky Way is the local escape speed v a . An estimate of v c 
may be obtained by studying the high velocity tail of the 
stellar velocity distribution in the Solar neighbourhood (e.g . 
iLeonard fc Tremainj|l99ol: IKochanekl 11999 : iMeiUod 11999) . 
Current samples include only a few tens of high veloc- 
ity st ars JCarnev et al.lll994l: lMeillonlll999l: ISakamoto et alJ 
120031) which only weakly constrain the shape of the velocity 
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Figure 3. Distributions, in real space and velocity space, of stars resulting from the accretion of a satellite galaxy by the Milky Way, 
illustrating the importance of radial velocities in the identification of stellar streams in the inner halo. Left: results for a stream in the 
outer halo (top panel: real space; bottom panel: velocity space) - the blue points represent stars within 10 kpc of the Sun; Middle: 
results for a stream in the inner halo (top panel: real space; bottom panel: velocity space); Right: distribution of integrals of motion for 
the inner halo satellite at the end of the simulation (top panel) and distribution of integrals of motion for stars originating in 33 individual 
satellites which could contribute to the stellar distribution of the inner halo (bottom panel). The data in each panel are convolved with 
expected Gaia observational errors in all quantities. See text for a detailed discussion. 



distribution. Radial velocities are necessary for this work as 
any underestimate of the space velocities of the stars trans- 
lates directly into an underestimate of the escape velocity. 
Further, given that the estimate of v e currently requires as- 
sumptions to be made about the shape of the velocity distri- 
bution of weakly bound stars, large samples of high velocity 
stars with accurate proper motions and radial velocities are 
essential to specify the shape without recourse to models. 
The accuracies required are not particularly demanding as 
the typical high velocity stars are moving at velocities in ex- 
cess of 350 km s -1 . BHB stars belonging to the stellar halo 
are a valuable tracer population for this work as they are 
numerous and relatively bright (Mv ~ 0.6). For a typical 
metal-poor, halo BHB star, the RVS will be able to obtain 
velocities to a limiting magnitude of about 16.5, correspond- 
ing to distances of about 15 kpc. There will be about 10 000 
such stars in the observable volume, which will place strong 
constraints on the actual shape of the local halo velocity 
distribution. 

The question of the full extent and mass of the Milky 
Way halo remains an open issue. Rece nt estimates (e.g. 
lKochaneklll996l;IWilkinson fc Evanslll999l) which include all 
available data on tracer objects outside 20 kpc still have un- 
comfortably large error bars. A radial velocity accuracy of 
10 — 15kms _1 is sufficient for this work. Stars at the tip 
of the RGB will be observable to distances of about 50 kpc 
while asymptotic giant branch (AGB) stars (Mv ~ —2.5) 



will be observable out to about 60 kpc. Samples of several 
hundred of each of these tracers will be observed by the RVS 
out to the orbit of the LMC. In addition, CH-type carbon 
stars with My ~ —2.5 will be observable out to 60 kpc. 
The density of these stars on t he sky is about four times 
higher than that of AGB stars llTotten fc Irwinlll998t) and 
so a sample of about 1000 objects is expected. 

Simulations show that a sample of 500 tracers extend- 
ing to large radii can reduce the error on the estimate of 
the enclosed mass to about 20 per cent and remove sy stem- 
atic errors (see Fig. 12 of IWilkinson fc Evans! Il99sl) . The 
RVS will provide such a sample extending to about 50 kpc. 
More importantly, it will also be possible to look for corre- 
lations in the motions of the tracers which would provide 
information about the origins of the stellar halo. Recent 
observations of stars in the outer halo suggest that a sig- 
nificant fraction of this Galactic component may be com- 
posed of streams jMaiewskil|2004l) . Knowledge of such cor- 
related motions is important for estimating the mass of the 
Galaxy as mass estimators generally assume that all trac- 
ers are drawn independently from an underlying velocity 
distribution. The unambiguous characterisation of a stream 
requires all three components of the stellar velocity - radial 
velocities alone place much weaker constrai nts, except for 
stars near the turning points of an orbit (e.g. IClewlev et alJ 
120051) . Gaia parallaxes will not provide accurate distances 
for these distant halo tracers. However, the Gaia-calibrated 
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distribution of absolute magnitudes for particular tracers, 
including any dependence of absolute magnitude on metal- 
licity or other internal parameters, will make it possible to 
obtain reliable photometric distances. If we conservatively 
assume that these distances can be calculated with an accu- 
racy of about 10 per cent, transverse velocities for all the 
above halo tracers can be obtained with uncertainties of 
about 10 — 15km s -1 . Thus, we can expect to have full space 
motions for samples of several thousand stellar tracers out 
to 50 — 60 kpc. Obtaining all-sky radial velocity coverage to 
an equivalent magnitude limit from a ground-based survey 
would be extremely time-consuming. 

The nature of the dark matter which makes up the 
dark haloes of galaxies is also an open question. In addi- 
tion to improving our knowledge of the distribution of mass 
within the Milky Way which itself constrains the proper- 
ties of the dark matter, the RVS will als o provide a direct 
test of the lumpiness of the dark halo. lYoo et all ||2004|) 
investigated the expected distribution of wide stellar bina- 
ries in haloes composed of massive compact objects (MA- 
CHOs) and found that encounters with MACHOs tend to 
disrupt the widest binaries and lead to a cut-off in the dis- 
tribution of binary separations. Using the observed distribu- 
tion of binary separations from a large sam ple of halo bina- 
ries JChaname fc Gouldll2004) . lYoo et aH conclude that the 
dark halo cannot contain a significant fraction of its mass in 
the form of compact objects with masses above 43M0 . The 
velocities furnished by the RVS, in combination with the 
Gaia proper motions, will permit confir mation of t he phys- 
ical association of all the binaries which lYoo et al] consider 
by providing the full space motions of their components. 
Further, it will be possible to investigate the Galactic orbits 
of these binaries which has implications for their expected 
survival times in a lumpy halo. In addition, a vastly larger 
sample of wide halo binaries will be obtained which will 
make it possible to apply this test of the nature of the dark 
matter over an increased volume of the dark halo. 

2.4 Dis-entangling the thin and thick disc 
populations 

The study of the kinematics of the disc of the Milky Way 
is important for a number of reasons. First, the thin and 
thick discs are major components of the Galaxy and pro- 
vide an opportunity to study the internal dynamics of a 
galactic disc at a level of accuracy which is impossible for 
external galaxies. Secondly, information about the forma- 
tion processes which produced the disc can be inferred from 
the motions of its stars. The RVS will be an essential tool in 
the study of the disc as it will provide velocities of sufficient 
accuracy to probe all aspects of disc kinematics. 

In order to investigate the formation of the Galactic 
thin and thick discs, it is vital to obtain precise information 
on the global Galactic kinematic properties (i.e. the velocity 
ellipsoids) so that deviations from expected behaviour can 
be reliably derived. At present, the kinematics of thin disc 
stars are poorly known: the structure of the velocity ellip- 
soid, its vertex deviation and inclination with respect to the 
Galactic Plane rest on a small sample of local stars. The ver- 
tex deviation measures the extent to which the local velocity 
ellipsoid is radially aligned and is strongly affected by the 
presence of moving groups in the local sample (see § 12. 5H . 



ISoubiran et alj i2003fl suggest vertex deviations which range 
from about 25 de grees for young stars to ab out degrees for 
older stars, while iDehnen fc Binnevl (Il998f) find that older 
stars have a vertex deviation of about 10 degrees. Addition- 
ally, to date, the vertical tilt parameter \(R) of the thin 
disc velocity ellipsoid remains ill-determined. The value of 
A indicates whether the velocity ellipsoids are aligned with 
the coordinate axes of a spherical (A = 1) or cylindrical 
(A = 0) coordinate system. In the first case the velocity 
ellipsoid points towards the Galactic center, while in the 
second it is everywhere parallel to the Galactic plane. The 
value of A is, in fact, strongly related to the coupling of the U 
and W velocities 4 and can only be properly constrained us- 
ing three-dimensional velocitie s. It is found to vary between 
0.4 and 0.7 at th e Solar circle dCudderford fc Amendtll99ll : 
lBienavmelll999l) . while no information is available at more 
distant locations. Although the value of A is directly related 
to the shape of the gravitational potential, the current range 
of values are consistent with a wide range of mass distribu- 
tions including an exponential disc mass distribution em- 
bedded in either a spherical halo with a flat rotation curve 
(A = 0.7) or a prolate halo (A < 0.5). More precise deter- 
minations of A, including its variation with position in the 
disc, are required and will be provided by Gaia. 

In addition to the alignment of the velocity ellipsoid, 
the magnitudes of the components of the velocity disper- 
sion tensor may also v ary according to th e age o f the stellar 
population considered. iLewis fc Freeman! lll989l) argue that 
a gradient in velocity dispersion as a function of stellar age 
ought to exist. The age-velocity dispersion relationship is 
based on samples o f stars in the Solar vicinity. For example, 
iQuillen fc Garnetl] (1200 if) find that the heating of the disc 
saturates after about 2 -3 Gyr. However, recent results by 
iNordstr om et alJ J2004h suggest a continued heating of the 
disc. The question is still under debate and larger data sets 
are required in order to resolve this issue. 

Information about the formation processes which pro- 
duced the thin and thick discs can also be gleaned 
from the rotational character of the discs. The pres- 
ence, if any, of vertical velocity gradients in the stel- 
lar components might be considered as another relict of 
the fo rmation proces s dChiba fc Beersll2000t ISoubiran et all 
20031: iGilmore et alj f2002l: iBeers fc Sommer-Larsenl Il995l 
Freenran^^Bl^md^Hawthor^^ iV-body simulations 

have demonstrated that radial mixing (e.g. due to the pres- 
ence of transient spiral arms) can wash out chemical and 
age gradients close to the Galactic pla ne on time scales 
of a few Gyr JSellwood fc Binnevl Eooil . Thus, it is pos- 
sible that only weak gradients will be detectable in the thin 
disc population. However, the thick disc should preserve a 
memory of its formation process since radial mixing is not 
so effective outside the Galactic plane and vertical gradi- 
ents in metallicity and velocity dispersion can, in principle, 
trace the formation history. These gradients are currently 
very poorly constrained. On the one hand, velocity gradi- 
ents perpendicular to the Galactic Plane are excluded by 
ISoubiran et al.l (120031) : on the other hand, it has been sug- 



4 The U direction is defined to be positive in the direction of 
the Galactic centre, while the W direction is positive towards the 
North Galactic Pole. 
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gested that the velocity dispersion and rotational velocity 
both decrease far from th e disc, with a veloci ty gradient 
of about 30 kms -1 kpc" 1 ijChiba fc Beersll200(t . The RVS 
will make a substantial contribution to the resolution of this 
issue by providing velocities for a large sample of tracers 
covering a significant fraction of the low density regions of 
the thick disc. For example, the velocities of K1III giants 
(with Mv = 1) will be measured to an accuracy of about 
25kms~ 1 up to distances of 25 kpc from the Sun. Prelimi- 
nary simulations show that Gaia will be able to detect the 
presence of a gradient in the thick disc rotation velocity of 
the order of lOkms -1 kpc -1 feertelli et al.ll200^ . 

It is also likely that the disc components of our 
Galaxy contain debris from past minor mergers. In par- 
ticular, the thick disc could be the product of the 
heating of a n ancient thin d i sc by a relatively large 
satellite (e.g. iQuinn et alJ 11991 [Velazquez fc White! Il99a 
iFreeman fc Bland-Hawthornl 120021) . Another related possi- 
bility is that the thick disc is constituted only by debris 
from disrupted satellite s having initially low inclination or- 
bits (lAbadi et alj 120031) : simulations show that a substan- 
tial fraction of the old disc stars could have formed in ex- 
ternal systems. Finally, a multi-component structure with 
tra cers of both mergers and accretions has been suggested 
by iGilmore et afl i2002fl . In all cases, kinematic substruc- 
tures are expected to remain in the disc and these should 
be detectable by the RVS, as well as in large radial velocity 
surveys such as the RAdial Velocity Experiment (RAVE: 
ISteinmetzl l2003fl . although the latter will survey a much 
small fraction of the disc than that covered by the RVS. 

Perhaps the earliest example of such an old kine- 
matic structu re with thick disc kinematics is the Arc- 
turus group jEggenl I1971T). whose nature needs to be 
confirmed ([N avarro et al 1 1200411 . The recent survey of 
IGilmore et alj (j2^0^i ~probed more distant regions of the 
thick disc, and found other hints of substructure. However, 
without full phase-space coverage for stars over a large frac- 
tion of the disc it remains difficult to make sense of the dif- 
ferent moving groups, and to establish whether or not they 
have a common progenitor. 

The outer disc of the Galaxy is very complicated, 
with observational evidence both for the warping and 
flaring of the stellar distribution outside the So lar cir- 
cle (e.g. lDiorgovski fc Sosinlll989l: lEvans et al.lll99Sl) . Kine- 
matic studies to date pr esent a confusing picture of the dy- 
nami cs of the warp (e.g. iMivamoto et al.lll988t ISmart et alj 
Il99ct) . The origin of the warp is also unclear with sug- 
gested possibilities including interactions between the disc 
and either the hal o, an infalling satellite or infalling gas (see 
iRobin et all 12001 for a discussion). While stellar proper 
motions are most suited for studies of warp kinematics, the 
RVS radial velocities can d etect the bulk motion associ ated 
with the warp. For example. lFeitzinger fc Spickerl (I1987T) dis- 
cuss the foreshortening effect of the warp on the distribu- 
tion of gas velocities in the Galaxy. For a star at £ = 90° 
(the direction of the maximum height of the warp above the 
Galactic plane) moving on a circular orbit, the foreshorten- 
ing effect amounts to a change in the radial velocity of the 
star of about 4kms _1 . The line of sight component of the 
rotation about the nodal line of the w arp is also comparable 
in magnitude (IMivamoto et al1ll99l) . 

Fig. 21 presents the results of a Monte Carlo simula- 
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Figure 4. Distribution of velocity residuals at various line-of- 
sight distances in the direction (£, 6) = (270, 45) when Gaia accu- 
racies on proper motions, radial velocities and distances are taken 
into account. The distributions were calculated from a Monte 
Carlo simulation of about 10 000 stars of thin disc, thick disc and 
halo populations brighter than V = 17. The three components of 
the observed velocity (radial, v^; vertical, ri z ; azimuthal, v^) are 
shown separately. 
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tion of about 10 000 stars of thin-disc, thick disc and halo 
populations brighter than V = 17 at ( L b) = (270, 45) de- 
rived using the Padova Galaxy M odel jBertelli et alJll99.4 
IVallenari et all2003t iBertelli et all2005t including expected 
Gaia accuracies on proper motions, parallaxes and radial 
velocities 5 . The differences between the observed and true 
velocity are plotted. The errors increase with distance, be- 
ing less than lkms -1 closer than 1 kpc, lOkms -1 at 4 kpc, 
and 20 — 30 kms -1 at 10 kpc. This accuracy should be com- 
pared with the typical velocity dispersion of the thin disc 
(less than lOkrns -1 ), thick disc (40 — 70 kms -1 ) and the 
halo (100 — 200 kms -1 ), respectively. Thus, it is clear that 
the velocity ellipsoids of both the thin and thick discs, as 
well as all other kinematic characteristics of their popula- 
tions, can be derived with unprecedented accuracy. 

The spiral arms associated with spiral density waves in- 
duce systematic motions of stars and gas both along and 
across the arms themselves and in the inter-arm regions 
which can be u sed to map out the gravitational poten- 
tial of the arm <Lin et al]|l96Sl) . Observations of OB as- 
sociations in the Cygnus-Orion arm have shown that the 
magnitude of these systemati c veloc iti es is between 10 and 
20 km s" 1 JSitnik fc Mel'Nikl Il999l) . iFeitzinger fc Spickerl 
il987l) found similar velocities in their study of simulated 
gas velocity fields in the disc. Thus, individual radial veloc- 
ities with accuracies of 5 kms -1 are sufficient to map out 
the velocity field associated with a spiral arm. The RVS 
will be able to obtain sufficiently accurate velocities of B5V 
stars to distances of about 2.5 kpc and of Cepheid tracers 
to distances of 6 — 10 kpc. These data will facilitate detailed 
comparisons with the spiral arms seen in numerical simu- 
lations and may provide clues to t he true origin o f spiral 
structure which is currently unclear iSellwoodll200bT) . 



2.5 Moving groups 

It is now well established that the distribution of stellar 
velocities in the Solar neighbourhood is not smooth - well- 
defined groupings in velocity space are clearly visible, for 
example in the proper motion d ata from Hipparcos (e.g. 
iFamaev et al.l2005t lDehnerll99Sf) . The clumping is not ran- 
dom as similar features are seen in the velocity distribu- 
tions of stars of different colours. The similarity between 
the colour-magnitude diagrams (CMDs) derived from stars 
which are members of moving groups and those of open clus- 
ters led to the suggestion that these groups were the rem- 
nants of dissolved clusters whose motion remains correlated 
long after the spatial coherence of the cluster has disap- 
peared. However, the wide age-range observed for stars in 
certain moving groups is difficult to understand in this pic- 
ture, suggesting that at least som e groups must have had 
a different origin. Recent work bv lDe Simone et alJ l|2004h 
suggests that certain old moving groups may be the product 
of disc heating by temporal variations in the disc potential 

5 Gaia is expected to measure parallaxes with a precision /n < 
10 per cent for a bright star having My = up to distances of 
10 kpc and with a precision of 1 per cent up to 2.7 kpc. Proper 
motions will be known with a precision of a few /ias/yr for stars 
brighter than V = 14 — 15 and of about 30/xas/yr for stars of 
V = 18, depending on the spectral type. 



due to the passage of stochastic spiral density waves. This 
would explain why stars of a variety of ages may be found 
in a single velocity structure. Outward moving streams (e.g. 
the Hercules stream) may be the result of chaotic relaxation 
of stars in the non- regular pa rts of phase space created by 
the central bar fe.g. lFuxll200f[ ). In addition, the outer Lind- 
blad resonance of the bar is near the Solar radius, w hich may 
give rise to some of the observed velocity features (|DehnerJ 

The characterisation of moving groups is of great inter- 
est as it provides clues to many aspects of the internal struc- 
ture of the Milky Way disc. The presence of these groups 
complicates the determination of the global properties of 
the disc populations. For example, the inclusion of a moving 
group in the local velocity sample can mimic the presence 
of vertex deviation of the velocity ellipsoid. Even the deter- 
mination of the solar motion relative to the local standard 
of rest has been shown to depend sensitively on whether or 
not all the nearby moving groups are i ncluded in the stella r 
sample used to measure this motion dFamaev et alJl2005l) . 
As was discussed above, the complete characterisation of the 
velocity distribution of disc stars is an essential step towards 
the development of a full understanding of the origin and 
evolution of the Milky Way disc. It will be greatly facilitated 
by the availability of radial velocities: the determination of 
the local velocity distribution using proper motions alone is 
equiv alent to performing a difficult de-projection jPehnenl 
Il998l) . The addition of the RVS radial velocities will make 
it possible to determine the full (i.e. unprojected) velocity 
distribution function for known groups. The identification 
of new groups will be made more robust through the use of 
integrals of motion to define membership rather than merely 
clumping in velocity space. The velocity accuracy required is 
modest: individual velocity errors of about 5 kms -1 will suf- 
fice to distinguish between typical moving groups. As Fig. 2] 
shows, moving groups will be easily identified out to dis- 
tances of around 2.5 kpc. 

2.6 OB Associations 

In recent years, observations of star forming regions have 
shown that stars almost never form in isol ation - most, if no t 
all, stars form in groups or clusters (e.g. lLada et al.lll99ll) . 
Thus, a full understanding of the details of star formation 
requires robust observations of young stellar clusters. The 
Milky Way disc contains a large number of clusters on a 
range of spatial scales from the low density, unbound and 
expanding OB associations (e.g. Perseus OB2) to the mas- 
sive open clusters (e.g. the Hyades and Pleiades clusters). 
The Gaia mission, and the RVS in particular, will improve 
our understanding of these objects by allowing the accurate 
determination of association/cluster membership based on 
distance and three-dimensional kinematics. The RVS will 
also facilitate measurements of the internal motions of stars 
within these systems. In addition, their collective space mo- 
tions within the Galactic potential will yield insights into 
the connection between star clusters and the main stellar 
populations of the disc. 

OB associations are unbound collections of recently 
formed stars which allow us to probe the initial mass 
function, primordial binary fraction, and early-time dy- 
namics resulting from the process of star formation. They 
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are ideal locations in which to determine the fraction of 
stars which form in binaries as they are young stellar sys- 
tems in which dynamical effects have had relativel y little 
impac t on the primordial binary population (e.g. iBrownl 
l2001bl) . Their internal velocity dispersions are typically a 
few km s~ . Membership has traditionally been determined 
using the convergent point method (which uses only po- 
sitions and proper motions but not parallaxes) although 
recent s tudies have intro duced modifications to this ap- 
proach Jde Bruiinelll999al) as well as new techniques such 
as the "spaghetti method" ijHoogerwerf fc Aguilarl Il999t 
lAguilar fc Hooeerwc ri 2001) which take account of parallax 
information. 

It has long been known that the intrinsic expansion 
of OB associations resulting from gas loss cannot be de - 
termined using proper motion data alone feraauwlll964l) . 
In particular, it is impossible to distinguish between a ra- 
dially expanding group at rest with respect to the ob- 
ser ver and a cluster app roaching along the line of sight. 
As lPerrvman et al.l l|l99Sfl discuss, the standard convergent 
point method for the identification of members relies on as- 
sumptions (small internal velocity dispersion, no internal 
kinematic structure or rotation) regarding the space motion 
of the cluster which can only be tested once radial velocities 
are available. For example, a net rotation may be inferred for 
a cluster with a large internal velocity dispersion if proper 
motions alone are used to determine membership due to the 
fact that the proper motion selection will produce an artifi- 
cial flattening of the cluster in velocity space. The presence 
or absence of rotation in star forming regions may provide 
important clues to the origin of the angular momentum ob- 
served in some globular c l usters (e.g.lMev lan fc Heggieil997t 
Ivan Leeuwen et alj|2000l : lAnderson^^Gn^]200^i and it is 
therefore important to be able to determine whether or not 
apparent rotation is real. The accurate determination of the 
orbits of young clusters and associations is of importance for 
the investigation of their relation to the surrounding field 
star population. In this context, the decomposition of any 
apparent radial velocity into bulk motion and expansion is 
clearly important. Only radial velocities, with accuracies of 
a few kms -1 , can resolve the situation. 

Radial velocities are an essential complement to proper 
motions in the correct identification of members of OB as- 
sociations which is the first step towards their use in discus- 
sions of the issues raised above. For example, using radia l 
velocities with errors ~ 3 km s" 1 . ISteenbrugge et al.l (12003ft 
have shown that in the case of the association Perseus OB2, 
a number of interlopers we re included in the list of possi- 
ble members determined bv lde Zeeuw et all Jl999Tl based on 
Hipparcos parallaxes and proper motions alone. The inclu- 
sion of radial velocity information in the analysis provided 
an additional selection criterion which made it possible to 
identify these interlopers as unrelated field stars, due to the 
lOkms -1 offset in radial velocity between the association 
and the disc stars along the line of sight. More accurate as- 
sessment of membership probabilities for stars near associa- 
tions is important as it allows uncontaminated CMDs to be 
plotted which in turn provide important information about 
the initial mass function. 

Associations are often studied using only their brighter, 
early type members (the OB stars). The intrinsic difficulty 
of obtaining radial velocities for early type stars, particu- 



larly when only the spectral region around the Call triplet 
is available, is well known. In the case of the RVS, for a 
slowly rotating B5V star (vsini = 50 kms -1 ) with absolute 
magnitude My = — 1, the expected velocity errors are about 
2.5 kms -1 at 1.5 kpc (V = 10) and about 5 kms -1 at 2.5 
kpc (V — 11), in the absence of reddening. Although radial 
velocities at this level of accuracy are insufficient to investi- 
gate the internal dynamics of associations (due to their small 
intrinsic velocity dispersions) , they will be useful for the de- 
termination of association membership. In addition, spectral 
information allows the estimation of the extinction towards 
individual stars which gives a further clue to whether or not 
a given star should be associated with a particular cluster. 

Associations also co ntain large nu mbers of lower mass 
(M < 2M Q ) stars (e.g. lBrowr]|2001bl) for which the RVS 
will easily furnish radial velocities at the 1 — 3 kms -1 accu- 
racy level. The exact numbers in particular associations will 
depend strongly on the local levels of extinction. However, 
we can estimate the range of masses which the RVS will 
be able to study using the properties of two known associ- 
ations. The nearest OB association is Scorpius OB2, which 
has a distance modulus of 5.3 — 5.8 magnitudes (equ ivalent 
to distance s of 116 — 1 44 pc; Ide Zeeuw et a lj|l_999|), solar 
meta llicity (lEggenll998T) . an age of 5-15 Mvr lde^3eus et alJ 
I1989T) a nd visual extinction in the range Av = 0.1 — 1.3 mag- 
nitudes <lde Bruiineil999bT) . The Cepheus OB3 association is 
one of the most distant associations currently known, with 
a distance modulus of 9.65 magnitudes (distance = 851pc), 
typical extinction Ay = 2.81 magnitudes, an ag e of less than 
10 Myr and solar metallicity iPozzo et all200af). We use the 
Padova stellar isochrones of iG^rardi^t^rtT feOOOF) assuming 
solar metallicity (Z = 0.02) and an age of lOMyr to deter- 
mine the stellar mass corresponding to the faintest magni- 
tude at which the RVS will obtain velocities to an accuracy 
of 3kms _1 . For K1III stars the magnitude limit is V = 16 
and for F2II stars it is V = 14. In the case of Scorpius 
OB2, we find that the RVS will be able to obtain velocities 
for stars with masses down to 0.4 — 0.7Mq depending on the 
extinction. In the more distant Cepheus OB3 this mass limit 
increases to about 2.2M© due to both its larger distance and 
higher levels of extinction. The availability of velocities for 
the lower mass stars in associations will allow more precise 
determination of the mean motion of the associations which 
in turn will permit better discrimination between members 
and non-members of all masses. 

The RVS data set will dramatically increase the sam- 
ple of OB associations which can be studied in detail. There 
are currently about 31 known OB associations within 1.5 
kpc of the Sun and which will be easily accessible to the 
RVS. Of these only 12 have been studied using a combi- 
natio n of Hipparcos data and ground-based radial veloci- 
ties ide Zeeuw et alJll999l). I n fact , 6 of the associations 
identified bv Ide^eeuw et al.l il999Tl were not included in 
previous catalogues which illustrates both the difficulty of 
identifying associations without full kinematic information 
and the likelihood that the number amenable to study with 
the RVS will be larger than expected. 

2.7 Open Clusters 

Somewhat more massive than OB associations, with grav- 
itating masses sufficient to retain a bound remnant follow- 
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ing the expulsion of the gas left over from star formation, 
the open clusters are important tracers of the young and 
intermediate-age stellar populations of the Galactic disc. 
Their value derives from the increased accuracy with which 
the distances, ages and metallicities of clusters can be de- 
termined compared to those of individual stars. Observa- 
tions of a large sample of open clusters will make it possible 
to look for metallicity gradients in the Galactic disc which 
may be preserved more strongly in the open cluster distri- 
bution than in the distribut i on of individu al stellar metal- 
licities (e.g. iFriel et~al"l 120021: rBrownll2001ah . Their internal 
velocity dispersions are generally a few kms -1 although 
some ( e.g. the Hyades clu s ter) c an be less than 1 kms~ . Re- 
cently, iKharchenko et aD (j200flT ) have used an all-sky stellar 
catalogue to determine the observed properties of a sample 
of 513 open clusters (each containing at least 18 members) 
within about 4pc of the Sun. The limiting magnitude in this 
study was V ~ 14 and the authors estimate that their clus- 
ter sample is complete to a distance of about 1 kpc. 

Gaia will be able to observe many open clusters - all 
known clusters will be amenable to study and it is expected 
that several thousand more will be identified within 5 kpc 
of the sun. The observability of an individual cluster de- 
pends on many factors such as age, distance and extinction, 
all of which affect both the numbers of stars with apparent 
magnitudes brighter than V = 17, and the velocity of the 
cluster relative to the local standard of rest, which can influ- 
ence the level of contamination by foreground stars. For ob- 
jects whose motion lies mostly along the line of sight, stellar 
proper motions will be small and radial velocities will place 
the tightest constraints on the membership lists. Thus, as 
in the case of OB associations, radial velocities are essential 
for the construction of an accurate census of members. 

For the Hyades cluster there are almost 400 stars 
brighter than V = 17 a t its present location of 
46.34pc iPerrvman et alJll998l) . For a similar cluster at a 
distance of 500 pc some 200 stars would still be observable 
by the RVS. At this distance the accuracy of the proper mo- 
tions will be about 0.5 kms -1 and so the goal of achieving a 
comparable radial velocity error with the RVS is achievable 
for only a tiny minority of member stars. For these stars, 
however, it will be possible to trace their orbits backwards 
in time and hence estimate a kinematic age for the cluster 
by determining their epoch of smallest separation. 

For nearby open clusters the RVS will provide ra- 
dial velocities with errors below lkms -1 for large num- 
bers of member stars. Even for the more distant clus- 
ters, however, the offset between the radial velocity of 
the cluster and the radial velocities of field stars along 
the line of sight means that the RVS radial velocities will 
generally be valuable in the construction of uncontami- 
nated membership lists. In addition, stars which have es- 
caped from the cluster but which remain on similar orbits 
within the Galaxy to that of the cluster centre of mass 
will be identifiable. These stars are vital to our under- 
standing of the disruption processes which affect open clus- 
ters (e.g. Ide La F uente_^laxros|[l99^JPOTte^ies Z wart et all 
l200ll) . For example. lde La Fuente Marcosl Jl996l) found that 
in TV-body simulations of open cluster dissolution the rem- 
nant of an open cluster which is left once the cluster has 
evaporated is very rich in binaries. Radial velocities will 
make it possible to confirm this result. The long term evo- 



lution of clusters in the Galactic disc is also of great interest 
as there is evidence that some of them, f or example Cygnu s 
OB2, are as massive as globular clusters llKnodlsederll200oT) . 
In this context, the existence of open clusters such as Berke- 
ley 17 with ages of about 9 Gyr dCarraro et a.1 j[T999l) allows 
us to investigate the long-term evolution of clusters in a 
strong tidal field. 

Young open clusters are valuable tracers of recent star 
formation and have been used to ide ntify the spiral arms of 
the Milky Way (e.g. lFeinsteinlll994T) . Their eventual disso- 
ciation builds up the population of field stars in the Galaxy 
and thus knowledge of the properties of as many such ob- 
jects as possible is an important step towards understanding 
how the stellar populations in the Galaxy formed and subse- 
quently evolved. As well as representing the dominant mode 
of current star formation in the Galactic disc, it has been 
suggested that the larger velocity dispersion of the thick disc 
compared to that of the thin disc might be partly due to the 
former containing a kinematically hot population of stars 
from clusters whic h became unbound by rapid gas expul- 
sion IIKroupall2002r i. The resultant characteristic mass func- 
tion of young clusters constitu tes an important test for this 
model faroupa fc Boilvll2002ft . A determination of the mass 
function of young clusters is also of fundamental importance 
if we want to understand the formation processes of clusters 
themselves - a complete census of cluster members for a 
large sample of young clusters is essential for this work. The 
cluster mass function also has implications for the chemi- 
cal evolution of our Galaxy and others because simulations 
have shown that long-lived, massive clusters can significantly 
enha nce the rate of Type la supernovae dShara fc Hurlevl 
120021) . 

The RVS will be able to test models of cluster disruption 
by looking for the unbound populations of stars surrounding 
young clusters which are the signature of the effects of gas 
expulsion. These stars will be easily identifiable due to the 
similarity of their metallicities and space motions with those 
of their parent cluster. In addition, the RVS will be able to 
detect the streams in the thick and thin disc populations 
which are expected to result from recently disrupted clus- 
ters. When a cluster becomes unbound, differential rotation 
and disc heating mechanisms spread the cluster stars within 
an elongated volume centred on the original cluster orbit. 
The total velocity dispersion of the stars increases with time 
- however, stars which are physically close together will 
have a lower velocity dispersion than t he initial value as a 
conse quence of Liouville's theorem (e.g. lBinnev fc Tremainel 
Il987l) . As Eq. shows, the velocity accuracy required to 
distinguish more than a few tens of streams is difficult to 
attain for a disc population with the dispersion of the thin 
disc. However, for an object which was disrupted during the 
past few hundred Myr there will be sufficient stars moving 
on orbits similar to that of the original cluster, and sharing 
common characteristics such as metallicity, to permit iden- 
tification of the remnant. 

2.8 Runaway stars 

Runaway stars are isolated, early type stars (spectral types 
O and B) with large peculiar velocities relative to the mean 
Galactic rotation. They are thought to originate in associa- 
tions and the feasibility of tracing their orbits backwards in 
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time in order to determine the cluster in which they formed 
has been demonstrated bv lHooeerwerf et al l fcOQcl) using a 
combination of Hipparcos data and radio observations. Two 
mechanisms which lead to the ejection of runaway stars have 
been suggested. One possibility is that close dynamical en- 
counters within an association can result in stars achiev- 
ing escape velocity. An alternative channel is the binary- 
supernova scenario in which the runaway star was originally 
a member of a close binary. Following the explosion of its 
companion as a supernova, and possibly the consequent dis- 
ruption of the binary, the runaway star acquires a sufficient 
velocity to move out of the cluster. Both mechanisms appear 
to occur in nature, but their relative importance ha s not 
yet been clearly established iHoogerwerf et all |200(J) . Re- 
tracing the orbits of large numbers of runaways is essential 
for the identification of the dominant production channel. 
This will, in turn, have implications for our understanding 
both of the internal dynamics of associations and clusters 
and of the fraction of high mass stars which form in bina- 
ries llPortegies ZwartH2000t) . ft will also improve our knowl- 
edge of the distribution of kick velocities acquired by pulsars 
at formation. 

Radial velocities are useful in this work for two reasons. 
First, knowledge of the full-space motion of the runaways 
makes the determination of their point of origin significantly 
more robust. Secondly, at present the identification of run- 
away stars relies on their large proper motions. The inclusion 
of radial velocity information will allow their identification 
based on a true space motion which may not lie in the plane 
of the sky - the RVS will therefore lead to a significant in- 
crease in the detectable sample of runaway stars. 

Another interesting possibility for the post- Gaia mod- 
elling of the Milky Way will be to use runaway objects as 
tracers of the gravitational potential of t he disc. An exam- 
ple of such an object is Cygnus-X2 which iKolb et aJl <l2000f) 
have suggested may have originated as an intermediate-mass 
X-ray binary in the Galactic disc which was subsequently 
driven to its present location 2.28 kpc out of the Galactic 
plane by the velocit y kick generated in a supernova explo- 
sion. More recently, iBrown et alJ <l2005l) have identified a 
hyper-velocity star moving at more than 700 km s -1 away 
from the Galaxy and currently located at about 55 kpc from 
the Galactic centre - its properties (including its age and 
metallicity) are consistent with its having been ejected from 
the Galactic bulge about 80 Myr ago. The Gaia data set 
will permit the identification of stars with anomalously high 
velocities compared to other stars in their vicinity. Deter- 
mination of the orbits of a significant sample of such stars 
using both radial velocities and proper motions may be used 
to place constraints (albeit somewhat circumstantial ones) 
on the distribution of mass in the Galaxy. 



2.9 Globular Clusters 

The globular clusters which orbit the Milky Way are valu- 
able laboratories in which to study the interplay between 
stellar evolution and stellar dynamics. Although the cen- 
tral regions of many Milky Way globular clusters will not 
be observable using the RVS due to crowding, the RVS will 
nevertheless contribute to our understanding of these sys- 
tems in a number of ways. For example, the outer regions of 



many clusters will be amenable to study and therefore the 
internal dynamics at intermediate radii can be investigated. 

It is reasonable to assume that the dense regions of glob- 
ular clusters will have a similar impact on the performance 
of the RVS to that of crowded Galactic fields of compara- 
ble surface brightness. As we discussed in Section ^. f I radial 
velocity estimates start to be degraded for stellar densities 
above 2 x 10 4 stars per square degree. For field stars, a limit- 
ing magnitude of V= 17 (i.e. a stellar density of about 2 x 10 4 
stars per square degree brighter than V = 17) corresponds 
to an integrated surface brightness (for all stars in the mag- 
nitude range V = 10 — 20) of V = 21.75 mag arcsec -2 (see 
IZwitteill2003allrl . Thus, stars with V= 17 will be observable 
in those regions of a globular cluster where the integrated 
surface brightness is fainter than 21.75 mag arcsec -2 . 

In order to determine which Galactic clusters will be 
observable, we assume further that a star with V< 17 will 
be observable against a given background if the difference 
between the magnitude of the star and the integrated sur- 
face brightness is the same as that between a V= 17 star and 
the V= 21.75 mag arcsec -2 critical surface brightness. This 
implies that, for example, stars in the ranges V= 14—15 and 
V= 15—16 can be observed against backgrounds of V= 19.75 
and V= 20.75 mag arcsec -2 , respectively. This assumption 
of a constant offset is probably conservative, since the S/N of 
the source spectrum increases for brighter sources - it may 
therefore be possible to measure velocities for the brighter 
stars against higher background d ensities. From the globular 
cluster catalogue of lHarria 1119961) we extract the concentra- 
tion and central surface brightness for each cluster and com- 
pute the radii within which the surface brightness is higher 
than the RVS limits for stars in a number of magnitude bins, 
assuming that the cluster can be represented by a King pro- 
file. We then calculate the proportions of the total number 
of stars in the cluster which lie outside these radii. Finally, 
using the distance, reddening, luminosity function and total 
mass of each cluster, we compute the number of observable 
stars in each bin. The cluster masses are estimated from 
their total luminosities using a linear interpolation between 
the 55 clusters for which mass estimates are available in the 
literature. The luminosity function of M92 is assumed for 
all the clusters. It is important to note that we have ignored 
the fact that for some RVS transits, individual spectra may 
overlap denser regions of the cluster, thereby reducing the 
effective number of transits per star. In view of this and our 
other simplifying assumptions, the estimates are probably 
only reliable to within about a factor 2-3. Nevertheless they 
give an indication of the likely contribution of the RVS to 
globular cluster studies. 

Fig. summarises the performance of the RVS in the 
Milky Way globular cluster population. In approximately 
33 clusters, the RVS will be able to observe more than one 
hundred stars brighter than V=15 to within 5 core radii - for 
24 clusters it will be possible to observe stars inside one core 
radius. The individual radial velocity errors for these stars 
will, in most cases, be smaller than the internal dispersion of 
the cluster, making them extremely useful for the discussion 
of the internal dynamics of the cluster. Thus the RVS will 
provide a large sample of stars for the detailed study of the 
internal kinematics of about 20-30 of the Galactic globular 
clusters. For 96 clusters, more than one hundred stars in the 
range V=15-17 will be observable to radii of less than 5 core 
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Figure 5. Performance of the RVS in Galactic globular clusters. 
For each cluster, the plot shows the innermost radius r (in units of 
the cluster core radius r c ) at which observations will be possible 
for stars of V= 15 (98 clusters: top panel) and V= 17 (121 clus- 
ters: bottom panel) and the total number of stars TV which will be 
observed in the cluster with V< 15 (top panel) and V= 15 — 17 
(bottom panel). Clusters which either cannot be observed or in 
which the innermost radius is more than 30 core radii have been 
omitted. See text for a detailed discussion. 

radii. These latter stars will be useful for determining the 
bulk properties of the clusters such as rotation. 

One key quantity which the RVS will constrain is the 
stellar binary fraction in globular clusters which, despite its 
importance for understanding the dynamical evol ution of 
clusters (e.g. lHut et al .111992 : 1 Wilkinson et aD^ro j j) . is cur- 
rently poorly constrained observationally (e . g lElson et alJ 
ll99alAlbrow et aljliooil: IStetson et al.ll2003l) . While we ex- 
pect that the fraction of tight stellar binaries should decrease 
with increasing radius in most clusters due to the effects of 
mass segregation, a limit on the binary fraction in the outer 
parts will at least provide a constraint on the fraction in 
the cluster as a whole. Cluster binaries will also be detected 
using the astrometric and photometric data from Gaia. 

Accurate velocities for stars in the outer parts of clus- 
ters will permit the study of rotation in clusters as well as 
providing estimates for the bulk radial velocities of the clus- 
ters. In conjunction with the Gaia proper motions for the 
clusters, these will determine the orbits of all the Milky Way 
clusters a nd place tighter constraint s on the mass of the dark 
halo (e.g. lwilkinson fc Evanslll999l. and §|2~31. 

The gravitational field of the Milky Way continually 
removes stars from the outer parts of the star clusters 
which orbit the Galaxy, leading to the development of tidal 
"tails" of stripped stars which extend along the orbits of 
the clusters. There is some evidence for th e presence of tails 
around a number of glob ular clusters (e.g. lLeon et al-lboOOt 
lOdenkirchen et al1 l2001). However, recent work has called 
into question the reality of certain of these tidal tails and, 
in particular, suggests that the tails of the massive globu- 
lar cluster Omega Centauri are artifacts produ ced by dif- 
feren tial reddening across the face of the cluster dLaw et alJ 
120031) . Kinematic studies of the outer parts of clusters are 
necessary in order to improve our understanding of cluster 
disruption by external tides. 

The RVS spectra will contribute to the resolution of this 
issue in two ways. First, as we discuss in § 12. Ill the spectra 
will facilitate the construction of accurate extinction maps 
throughout the Galaxy. Secondly, the velocities of stars in 



the tidal tails of a star cluster are similar to the bulk motion 
of the cluster and, in general, differ significantly from those 
of surrounding field stars. As a result, the velocity estimates 
provided by the RVS (and, of course, by the Gaia proper 
motions) will often be sufficient to distinguish between es- 
caping cluster members and field stars. Thirdly, the metal- 
licities of cluster stars will generally differ from those of the 
surrounding field stars. For brighter stars, this will provide 
another potential confirmation of a physical connection, if 
any exists, between the "tails" and the cluster. 

All Galactic clusters with sufficient numbers of bright 
(V< 17.5) stars will be amenable to study via a combina- 
tion of the above approaches; for example, the well-studied 
cluster Pal 5, whose tidal tails have now been traced in the 
stellar number density dist ribution to distances of about 2 
kpc from the cluster centre llOdenkirchen et al.l2003|) . While 
the internal velocity dispersion of the tails is too small for 
the RVS velocities to contribute to discussions of their in- 
ternal str ucture (the intrinsic dispe rsion of Pal 5 is less than 
lkms -1 llOdenkirchen et al.ll2002T) and the spatial coher- 
ence of the streams strongly suggests that they must also be 
kinematically cold), the velocities will be of sufficient accu- 
racy to distinguish tail members from Galactic field stars. 
This will be of particular value for those clusters whose tidal 
tails are less extended or less clearly defined spatially than 
those of Pal 5. It is also worth noting that Gaia will also 
contribute to the identification of tidal streams by means 
of proper motions, since stars in tidal tails display proper 
motions very similar to those of cluster members and, in 
general, q uite different from the field stars along the line of 
sight (e.g. lKing et"ai1ll99Sl) . 

2.10 Chemical Evolution 

To build a realistic model of the chemical evolution of the 
Gala xy one needs as many observational constraints as pos- 
sible jPageJll997l) . These observations not only have to in- 
clude the relative numbers of stars in the thin and thick 
disc or the halo, but also have to be able, using the stellar 
abundances of some chemical elements, to constrain partic- 
ular model inputs, for example the rates of type I and II 
supernovae. Possible by-products of these chemical abun- 
dance measurements could include the determination of an 
Age-Metallicity relation and possible gradients in the thin 
and thick discs (see also § 12.131 . The RVS spectra will re- 
veal mainly [a element/Fe] ratios for late type stars (e.g. 
the a elements Ca, Mg, Si, Ti), which constrain the sites of 
nucleosynthesis - these stars will require follow-up observa- 
tions, for example, of the r process elements not obtainable 
from the RVS spectra. The RVS will provide precise chemi- 
cal abundances for 3 — 6 x 10 6 stars brighter than V=12-13. 
This will make it possible to probe the thin/thick disc and 
halo populations up to a distance of 2 kpc when K1III stars 
having absolute magnitude Mv ~ 1 are used as tracers and 
up to 300 pc using G2V stars of Mv = 5. 

The variation of [a/Fe] among stars as a function 
of their location in the Galaxy will most strongly con- 
strain models of the quantity of gas in the Galaxy at dif- 
ferent epochs and also the supernova rate. The combina- 
tion of chemical abundances and kinematical properties 
has revealed some difficulties with nai ve infall models of 
the Galaxy as was first demonstrated bv lNissen fc Schuster! 
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Il997|) . It now seems difficult to support a simple infall 
model: this adds weight to the introduction of a more violent 
scena rio of the formation of the Galaxy fe.g lHelmi fc White! 
Il999l) . However, in such a scenario it is not obvious that 
one would expect a well-d efined relation betwe en chemical 
abundance ratios and age l|Rvan fc SmitblfeoOSF ) during the 
epoch when Type la supernovae are supposed to explode 
([Fe/H] ~ -1.0). 

Searching for extreme Population II stars is also a chal- 
lenge for the Gaia survey in order to determine the prop- 
erties of the most metal-poor stars. These fossil stars are 
representative of the first phase of the Galaxy just after 
the hypothetical era of the Population III stars. Observing 
their [C/O] ratios with ground-based spectroscopy makes 
it possible to look for evidence of the presence of Popula- 
tion III stars very ea rly on in the history of the Galaxy. 
lAkerman et al.l (I2004T) have discovered a possible increase 
of the [C/O] ratios for very metal poor stars which would 
require the presence of massive Population III stars at the 
beginning of the formation of the Galaxy. Calcium lines will 
be visible even in the RVS spectra of very metal poor stars. 
The Ca line at A = 854.2 nm can be seen with an equivalent 
width EW = 0.0482 nm at [Ca/H]= -3 and EW = 0.0245 
nm at [Ca/H]= -4 for a giant star with T cff = 4800 K, 
logg = 1.5. This should allow the separation of a giant with 
[Ca/H] = —4 from a giant of [Ca/H]= —3 down to about 
a S/N~ 10 corresponding to V~ 15 for the above types of 
star. Assuming this star is a giant of absolute magnitude 
My = — 1 then a survey for stars more metal poor than 
[Ca/H]= —3.5 will be possible within a sphere of 5-10 kpc. 
Gaia will also be able to detect binarity among such halo 
stars and hence rule out mass transfer scenarios through 
which the carbon abundance in the primordial atmosphere of 
an extreme-Population II star is modified, thereby produc- 
ing an incorrect diagnostic on th e presence of Populati on III 
stars in the early phase (see e.g. iMasseron et alJl2003l) . 



Among the successes of the standard big-bang model 
and of the W ilkinson Microwave A nisotropy Probe (WMAP) 
experiment iSpergel et alj 120031) are the stringent con- 
straints on predictions of primordial light element abun- 
dances. In particular, the helium abundance Y is now con- 
strained to lie in the range 0.247 — 0.252. This value is then 
supposed to be the primordial helium content of the Galaxy. 
Its enrichment with time is assumed to be a consequence of 
multiple generations of stars and is expected to follow the 
enrichment of heavy elements Z. However, the details of the 
helium to metal enrichment relation SY/SZ which describes 
the increase in helium from p rimordial levels have long been 
debated in the literature (see lHog et al.lll993 . for a review). 
Gaia will provide a survey of several thousand binaries many 
of which will have accurate age estimates because when the 
masses of the components of a binary are known with accu- 
rate effective temperature, luminosity and abundance Z then 
the age can be determined by fitting evolutionary tracks to 
the two error boxes of the binary in the HR diagram. The 
ages and estimates of He content deduced from these Hr di- 
agram fits, can be used to discuss the helium enri c hment 
of individual systems: for example, iLebreton et al.l (1200 J) 
present an analysis of five binaries in the Hyades cluster 
and determine the helium content and age of the cluster. 
The Galactic enrichment relation SY/SZ will be constrained 




5000 io 4 
leliocentric distance (pc) 



1.5x10 



Figure 6. Extincti on along several lines of sight towards the 
Galactic centre. The IVallenari et alj l2003f) determination (solid 
line in each panel) is the mean value in a field of 0.5° X 0.5°. The 
error bars indicate variable re ddening inside the field . The d otted 
line represents the model of iMendez (fe 1 _vaii_Altena[ _ 1 1998 ) : the 
dot-dashed line with triangles is the Erlnmiel^^SDe^eiril^OOlf) 
mo del including a r eseati ng factor. The star gives the value from 
the ISchlegel et all <1998Tl maps. When the determination from 
the ISchlegel et alJ maps lies outside the plotted range, a labelled 
arrow indicates the reddening value. 



not only for close binaries with Solar metallicities but also 
for very metal-poor binary stars. 

At the same time, Gaia will provide the opportunity to 
measure magnitudes, absolute distances and chemical com- 
positions for a substantial fraction of halo and thick disc low- 
mass Horizontal Branch (HB) and RGB stars. These new 
data will supply accurate estimates of the R paramete r, i.e. 
the r atio between the number of HB and RGB stars (|IberJ 
Il968t) . The comparison between empirical star counts and 
evolutionary lifetimes provides an estimate of the Helium 
content and in turn an upper lim it on the primordial Helium 
abundance iZoccali et alJl2000blh Up to now the R parame- 
ter has been estimated only in globular clusters and current 
estimates are partially hampered by sta tistics and by the ac - 
curacy of spectroscopic measurements JCassisi et alJl2003h . 
The new Helium abundances will also provide the opportu- 
nity to estimate SY/SZ in the Galactic halo and thick disc 
as a function of Galactocentric distance. 



2.11 Extinction maps 

Interstellar absorption represents one of the major compli- 
cations in the simulation of the CMDs of Galactic stellar 
populations along any line of sight, especially at low Galac- 
tic latitudes across the disc. The uncertainties in star count 
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predi ctions caused by r eddening may amount to some 26 per 
cent ilChen et aI.lll999Tl . 

Several empirical mode ls of Galactic absorpt i on are 
available in the literature. iMendez fc van Altenal il99cf) 
make use of the large-scale properties of the dust layer in 
the Galaxy to derive the absorption in the Galactic plane. 
ISchleeel et all (Il998h use the COBE /D1RBE 100 ^m and 
240 /im data to construct maps of the dust temperature 
in the Galaxy. In high-latitude regions, the dust map cor- 
relates well with maps of H I emission, but deviations 
are found in parts of the sky and are especially conspic- 
uous in regions of saturated H I emission towards denser 
clouds and in areas of H2 formati on in molecular clouds. Re - 
centlv. lDrimmel fc Soergel <200ll) and lDrimmel et alj (I2003T) 
presented a three-dimensional model of the dust distribu- 
tion based on COBE /DIRBE infrared data. As stated by 
iDrimmel et alj (l2003f) . regions having anomalous emission 
due to warm dust are not well described by the model. 

When comparing the various literature models of ex- 
tinction in analyses of CMDs observed along the line of sight, 
very large discrepancies soon become apparent. First, there 
are still problems with the accuracy of the zero point cali- 
bration of the extinction maps. In fact, as iBursteinl (l2003f) 
points o ut, even at high Galactic latitudes the reddening 
map by iBurstein fc Heilesl lll97ctl provides reddening val- 
ues i_that_are_0^2 inag smaller than the reddening map 
bv lSchleeel et all <ll998T) . At low Galactic latitudes and to- 
wards the spiral arms, the di fference may be as lar ge as 4-5 
magnitudes (see Fig. 151 and IVallenari et al.l I2003T) . This is 
due in part to variations of the dust properties on small 
scales and in part to the many uncertainties which are still 
associated with dust emission models. 

The RVS can directly measure an interstellar or circum- 
stellar extinction corresponding to Eb-v = 0.10 using the 
862 nm diffuse interstellar band (DIB: see Paper I) of early 
type stars brighter than V ~ 12 — 13 which has been shown 
to be directly associated with the dust phase of the interstel- 
lar medium. Simulated RVS data indicate (see Fig. [7| that 
the imprint of the 862 nm DIB will be detected even in the 
spectra of much fainter stars with magnitudes up to V ~ 16 
with sufficient accuracy to trace not only the distribution of 
the interstellar medium but also the radial component of its 
kinematic motion (the Doppler velocity of the mass center of 
the dust cloud is calculated from the wavelength position of 
the DIB line centre). The resulting three-dimensional maps 
of interstellar extinction are independent of the photomet- 
ric approach to reddening determination which is based on 
a comparison between observed colours and modelled intrin- 
sic colours. Thus, the RVS will make a vital contribution to 
the construction of a new generation of accurate star count 
maps of the Galaxy which are essential to the development of 
complete models of Galactic structure. In addition, accurate 
extinction-corrected stellar magnitudes will be of enormous 
importance for the testing of stellar models (see §0|l, solv- 
ing the degeneracy between reddening and temperature that 
usually affects determinations from photometry. 

As an illustration of the effects of reddening on 
the stellar distribution in different regions of the Galac- 
tic spheroid we selected a few well-observed fields. In 
particular, we focused our attention on t hree stan- 
dard stellar fields for which IStetsonl teOOCf) (see also 
: //cadcwww.hia.nrc . ca/standards/ I collected new 
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Figure 7. Recovery of the 862 nm DIB from simulated RVS 
spectra. The stellar spectra are o btained from the st ellar popula- 
tion synthesis Galactic model of iRobin et al 1 l2003f>. The three- 
dimen sional model of the dust distribution of lDrimmel fc Spergell 
l200lf> is used. Left: Number of simulated stars as a function of 
distance in three Galactic directions (I = 10°; b = 0°, b = 10° and 
b = 30°). Black represents all simulated stars, dark grey those for 
which the DIB line was recovered with an accuracy of at least 
30kms _1 and light grey those for which the recovery of the DIB 
line was more accurate than lkms -1 . As expected, extinction is 
a strong limiting factor only within the Galactic plane. Right: 
Number of simulated stars as a function of inaccuracy in equiv- 
alent width EW, i.e. the difference betwee n the simulate d and 
recovered EW of the DIB line. According to lMunaril <200Cl) . EW 
is proportional to Eb-v (an EW of 0.05 nm corresponds to a 
reddening Eb— v = 1.35). 



and homogeneous multi-band photometric data. Table[3]lists 
the coordinates and positional parameters of the selected 
fields. 

Fig. Ushows the (V - I, I) CMDs of the selected fields 
which are assumed to be representative of the number densi- 
ties and extinctions for intermediate latitude Galactic fields. 
The data plotted in this Figure clearly show the effects of 
reddening. When moving from bottom to top, the redden- 
ing increases from E(B-V)=0.04 to 0.16: the increase, as 
expected, causes a systematic shift towards fainter mag- 
nitudes and redder colours for the field stars. This means 
that observed luminosity functions and colour distributions 
are strongly affected by the interstellar extinction and by 
differential extinction, if any, along the line of sight. We 
note in passing the large di fference in the me a n inte rstel- 
lar extinction between the IBurstein fc Heilesl il982T) and 
the lSchlegel et ail dl99ctl reddening maps (see column 7 in 
Table[3J . This emphasises again the crucial role of individual 
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Field FOV (arcmin 2 ) RA(J2000) DEC(J2000) 1 b A B Ns 



L95 46.0 x 50.9 03 54 12.40 +00 13 09.0 188.79146 -38.25715 0.67/1.51 870/1090 

L107 44.1 x 37.5 15 39 20.90 -00 20 25.9 5.70550 41.22154 0.29/0.48 1170/1525 

MarkA 33.9 x 26.9 20 43 43.20 -10 46 16.9 35.89721 -29.77091 0.18/0.25 1700/2215 



Table 3. Parameters of selected Galactic fields which illustrate the range of conditions which the RVS will encounter at intermediate 
latitudes. The table presents the field name, the field of view for each observed field in square arcmin, the positions of the field centres in 
ecliptic and Galactic coordinates, the interstellar extinctions Ab estimated over an area of 5x5 arcmin centred on the individual fields 
and the estimated number of stars in each field N B which will be detecta ble by the RVS. Note: T wo v alues of Ar (and the associated 
N a estimates) are quoted for each field. The Ab values are taken from the lBurstein fc Heilesl Jl982T) and lSchlegel et all Il998f) reddening 
maps, respectively. 
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Figure 8. (V - I, I) CMDs for the three fields listed in TablcEl 
The dotted horizontal lines show the likely range of the faint 
magnitude cut-off for the RVS radial velocity determinations. 



reddening measurements in the improvement of the accuracy 
of stellar parameters. 



2.12 Star formation history 



composition and distance are directly exploitable. Addition- 
ally, the age-metallicity degeneracy further complicates the 
problem. Following the Hipparcos mission, it was possible 
for the first time to derive the CMD, in absolute magni- 
tudes , of field stars in the Solar vicinity iPerrvman et alJ 
I1995T and from this CMD to study the history of the Solar 
neighbourhood. 

The Gaia mission will significantly exceed the perfor- 
mances of Hipparcos permitting the determination of the 
star formation histor ies (SFH) of th e disc and halo of the 
Milky Way. Recently. iBertelli fc Nasi feOOll) determined the 
SFH of the Solar vicinity from Hipparcos data, covering the 
total lifetime of the disc (10 Gyr). They find that from the 
Hipparcos catalogue it is possible to select a complete sam- 
ple of stars with well measured parallaxes down to My = 4.5 
in order to include main-sequence and evolved stars, inside 
a sphere of radius r = 50 pc. Gaia will be able to observe a 
much larger sample of stars covering in distance and position 
a large portion of the Galaxy and possessing the same degree 
of accuracy as that obtained by Hipparcos. iBertellil (I2002TI 
has demonstrated that Gaia will indeed allow the study of a 
sample of stars complete down to My = 4.5 with a parallax 
accuracy better than a^/n < 0.1 up to a distance of 2 — 3 
kpc. Here it is worth remembering that due to the effects of 
radial mixing in the disc as a response to spiral wave per- 
turbations, old stars born at the location of the Sun may 
now be spread over a range of Gal actocentric radii from 4 to 
12 kpc JSellwood fc BinnevlEx)^ . Within this volume, the 
RVS will provide information about the average metal con- 
tent of the stars brighter than about V — 14—15 (see Tabled 
and Paper I), solving the age-metallicity degeneracy which 
always hampers the estimation of the star formation rate 
from Hertzsprung-Russell (HR) diagrams, in particular for 
older stars. Additionally, information on the metal content 
of fainter stars will be provided by Gaia photometry and 
astrometry. Since different populations have different kine- 
matics and metallicities, coupling these properties would al- 
low us to distinguish statistically between thin/thick disc 
and halo stars. 



The determination of the history of star formation from the 
CMDs of composite stellar populations is an important goal 
of modern astrophysics. The problem is easier to tackle in 
galaxies in which individual stars are resolved and CMDs 
are derived, because we may assume that all the stars lie 
at almost the same distance. However, in our own Galaxy 
the problem is significantly more complicated because there 
are differences in the distances of the stars and therefore 
only CMDs containing stars of inhomogeneous age, chemical 



2.13 Age-metallicity relation 

The age-metallicity relation (AMR) for disc stars, if 
such a relation exists, gives information about the pro- 
cess of stellar formation, stell ar orbit diffusion from scat- 
tering by molecular c louds iFrancois fc Matteuccil Il993l : 
lEdvardsson et alJ Il993ft and about the time scale for gas 
mixing lvan^CTiHoek&: de JonJl997l) . The existence of an 
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AMR has long been a controvers ial issue. iTwarod ( r98C ) 
found an AMR in field stars, while lEdvardsson et alJ ( 1995 ) 



derive no AMR from a sample of about 187 FG giants with 
known metallicity, distance and magnitude. iNg fc Bertellil 
lll998l) . using revised age estimates, derive a moderate AMR 
with a slope 0.07 dex/Gyr. These authors find that the main 
source of uncertainty in the age- metallicity determination is 
due to the distance estimates which are used in the conver- 
sion from apparent to absolute magnitudes. On the basis 
of their simulations, distances need to be known with at 
least 5 per cent accuracy to obtain ages with a precision 
of about 16 per cent. Gaia will observe stars brighter than 
Mv = 5 with a distance accuracy of less than 1 per cent up 
to 1 kpc, and with an accuracy of 5 per cent up to 2 kpc. 
The results of lNg fc Bertellil 1 19981 are i n substantial agree- 
ment with the lRocha-Pinto et al.N2000t) study of a sample 
of 552 stars. They derive metal content information from 
Stromgren photometry and ages from chromospheric activ- 
ity, finding an AMR of 0.05 dex/ Gyr. 

Based on 5800 stars from the Hipparcos catalogue with 
ages derived from isochron es and metallicities estimated 
from Stromgren photometry, iFeltzing et al.l (1200 if) found an 
AMR only for objects younger than 2 Gyr. Recently, using 
isochrone fitting to derive a ges and Stromgren pho tometry 
to estimate the metallicities. iNordstrom et al.l i2004T) find no 
AMR on a large sample of 14000 F and G dwarfs in the solar 
vicinity. All the above studies, however, suffer either from 
the small sample of well- measured stars fedvardsson et all 
ll99fllNg fc BerteUill998h or from uncertainties in the deter- 
mination of the metal content f rom Stromgren photometry 
which is difficult to calibrate (lHavwoodll2002l) . Gaia will 
allow age estimations via accurate positions in the HR di- 
agram. The determination of the metal content for a large 
sample of disc stars which will be possible using the spectra 
from the RVS will be an essential complement to these accu- 
rate ages in improving our knowledge of the age-metallicity 
relation. 

Finally, chromospheric activity can be calibrated as an 
age indicator against isochrone fitting. Chromospheric ac- 
tivity is expected to decline while stars are aging and can 
therefore be used to date stars, mainly F, G, and K dwarfs. 
At present, this method is poorly calibrated and many un- 
certainties remain: first, there is intrinsic variability of stellar 
activity (e.g. the activity cycle in the Sun); secondly, stel- 
lar activity is caused mainly by rotation which, although it 
generally decreases with age, can be infl uenced by, for exam- 
ple, tidal interaction in binary systems feawaleii | l989f) . The 
RVS data set will include thousands of eclipsing binaries (for 
which the uncertainty in the projection of the rotation vec- 
tor along the line of sight is removed) containing F,G and K 
dwarfs locked in synchronous rotation with their orbital mo- 
tion and spanning a wide range of orbital periods. It will thus 
be possible to calibrate accurately the correlation between 
their levels of chromospheric activity (as traced by photom- 
etry and Call emission cores in RVS spectra) an d their rota- 
tion speed (see § 14.311 . Using the lKawalerl il989T> calibration 
of rotational velocity against age for main-sequence stars 
older than 100 Myr, an uncertainty of 5kms _1 in the rota- 
tional velocity will result in an uncertainty of 20 per cent in 
the age of a star with a rotation velocity of 50kms~ . 
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Figure 9. CMDs for a field in Baade's window towards the glob- 
ular cluster NGC 6522 (top panel) and a region in the Blanco 
field (bottom panel). 



2.14 The Galactic Bulge 

The derivation of the ages, age distribution, metallicities 
and kinematics of the stellar populations of the bulge is ex- 
tremely important for our understanding of the process of 
galaxy formation. As was discussed earlier, the high stellar 
densities along lines of sight which pass within 5 — 10 de- 
grees of the Galactic centre render observations of the bulge 
using the RVS very difficult. However, the outer bulge will 
still be amenable to study by the RVS and much will be 
learned about the nature of the bulge/bar from the RVS 
observations. 

To obtain a more quantitative estimate of the number of 
stars appearing in the RVS band towards the central regions 
of the Galactic bulge we observed two regions, one located 
at I = 0.99, b = —3.94 (in Baade's window) which contains 
the Globular Cluster NGC 6522 and one at I = 0.27, b = 
—6.19 (the Blanco field) which is one of few regions of the 
Galactic bulge that is characterised by low reddening. We 
collected a series of CCD images in the Stromgren bands 
u,v, b,y, Ca with the Danish 1.5m telescope at ESO La Silla 
(field of view 13.7 x 13.7 arcmin 2 , 0.46 arcsec/pixel). Fig. EH 
shows the (b — y,y) CMD for these fields. In the Blanco field 
we detected 2700 stars brighter than y = 17.5, while the 
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Baade's Window field contains 8900 such stars (we note that 
the y band is the closest Stromgren filter to the RVS band) ; 
these correspond to densities of about 52 000 and 170 000 
stars per square degree, respectively. The mean reddening 
towards both these fields is E(B-V)« 0.45. The difference in 
the star counts is due to the presence in the latter field of 
NGC 6522 and also to the fact that in this region we are 
looking towards the Galactic center. 

As the Figure shows, the majority of bulge stars lie 
well below the likely magnitude limits of the RVS in these 
regions. Investigations are underway to determine strategies 
for obtaining at least some information on the brighter stars 
in certain bulge fields. In the event that it is possible to ob- 
tain these data, Fig.|5|shows that bulge AGB stars would be 
amenable to study at magnitudes between 14 and 15. Given 
that the line of sight depth of the bulge is approximately 
one magnitude, a limit of 14.5 would permit the observation 
of AGB stars all through the bulge on these lines of sight. 

Even in the absence of data close to the Galactic cen- 
tre, the RVS will illuminate certain aspects of the Galactic 
bulge, in particular the presence and nature of the Galac- 
tic bar. The presence of a bar may be the origin of some 
of the stellar moving groups observed in the Solar neigh- 
bourhood. Further, there is growing evidence for the pres- 
ence of bulge/bar features in the stellar distributions at 
larg er angular distance s from the Galactic centre. For exam- 
ple, |Sammera^^^^L| J2000T) identified a red clump feature 
at I = 27° and a similar one h as been identified at I = 9° 
bv lBabusiaux fc Gilmorel i2005l) . This feature, which i s also 
seen in observations of OH/IR stars JSevensterlll999D and 
SiO masers ilzumiura et al.lll999t) . may be associated with 
a secon d bar, such as are ofte n observed in external galax- 
ies (e.g. lErwin fc Sparkel20o1 . or a stellar ring at a distance 
of 3 to 4 kpc f rom the Galac t ic cen tre. Along the same lines 
is the work bv lPicaud et alJ ll2003f) who detect a star-count 
excess in NIR data with respect to the expected disc popula- 
tion at 15° < I < 21° . The radial velocities provided by the 
RVS will help to distinguish between the possible explana- 
tions of these features. Thus, although the RVS will at best 
observe a limited number of bright stars in the most central 
regions of the Galaxy (and at worst not observe the central 
regions at all), it will nevertheless provide very useful con- 
straints on the properties of the bar and bulge by observing 
bulge features outside 10° from the Galactic centre. 



2.15 Local Group Galaxies 

2.15.1 The Magellanic Clouds 

The Magellanic Clouds provide a unique opportunity to 
study the evolution of dwarf irregular galaxies at close quar- 
ters as well as the effects of galaxy interactions (both the 
interactions between the clou ds themselve s and the exter- 
nal effect of the Milky Way - ) . ISpitel J2002I) has highlighted 
a number of outstanding issues regarding the Magellanic 
Clouds which can be addressed by the Gaia mission. The 
resolution of one of these in particular relies on the spec- 
troscopic capabilities of the RVS: does the Large Magellanic 
Cloud (LMC) possess a pressure-supported stellar halo? Ev- 
ide nce for such a popu lation has recently been presented 
bv lMinniti et alJ i2003fl using observations of 43 RR Lyrae 
stars in the inner regions of the LMC. The observed velocity 



dispe rsion of this stellar halo is 50 ± 10km s _1 dMinniti et alJ 
2003), which is considerably larger than the LMC thick disc 
velocity dispersion of about 20km s _1 , estimated from the 
motions of the o ld LMC star clusters and intermediate- 
age Carbon stars (Ivan der Marel et al.ll2002l) . The RVS will 
therefore be able to distinguish kinematically between mem- 
bers of the two populations for the RR-Lyrae and metal- 
poor K giants with magnitudes in the range 16-17.5 as it 
will yield velocities for these stars accurate to better than 
15km s _1 . It will observe stars over the entire area of the 
LMC thereby making it possible to determine whether the 
population is rotating - this is not currently constrained by 
the data of lMinniti et al.H2003l) . The presence of such a halo 
has direct implications for our understanding of the forma- 
tion history of the Magellanic Clouds, as it implies a similar 
hierarchical evolution to that of the Milky Way. It also has 
implications for the interpretatio n of the mic rolensing data 
towards the clouds. For example, lAlvesI (120041) suggests that 
two of the known MACHO microlensing events might be due 
to lensing of LMC disc stars by stars in the putative halo, 
although he points out that the halo cannot account for all 
the observed microlensing optical depth. In addition, a com- 
parison of the kinematics of other stellar populations with 
those of the Carbon stars may yield further insights into 
the nature o f the stellar bar in th e LMC, which is currently 
unclear (e.g. IZhao fc Evansll200o|) . 

The Ma gellanic Clouds contain a large population of 
star clusters iMackev fc Gilmorei2003 'allpT) with integrated V 
magnitudes of about 10 — 13. As in the case of the Milky Way 
globular clusters, these are useful tracers of the mass distri- 
bution and internal kinematics of the Magellanic clouds. The 
RVS will be able to obtain accurate radial velocities for all 
the clusters yielding a data set of more than 50 tracers per 
galaxy. Given that the cluster populations extend to larger 
radii than other tracers, these will be invaluable in the deter- 
mination of the total gravitating mass of the clouds thereby 
constraining models of their formation. The internal dynam- 
ics of the Clouds and their relationship to each other, as well 
as the role played by the Milky Way in their evolution, are 
issues which remain to be comprehensively addressed. 

2.15.2 The Andromeda Galaxy 

The Andromeda galaxy is the other massive galaxy in the 
Local Group. Its stellar halo is currently the focus of consid- 
erable interest due to the presence of la rge amounts of sub - 
structure, in particular a stellar stream jlbata et alj|2001ah . 
Recent estimates of the total mass of the Andromeda galaxy 
have shown that there is no kinematic evidence for the gener- 
ally held belief that Andromeda is the most massive galaxy 
in the Local Group: the data on tracers outside about 20 
kpc favour a halo which falls o ff more rapidly outside 3 
kpc than that of the Milky Way fevans fc Wilkinso rfeOOOT) . 
The key difficulty facing attempts to verify this conclusion is 
the paucity of tracer objects in the crucial radius range from 
30 to 100 kpc. At present there are only a handful of globular 
clusters and dwarf galaxies at these radii. Recent observa- 
tions of the globular cluster population of M31 within 25 
kpc have shown that rough ly half of the metal- poor clusters 
are brighter than V = 17 JPerrett et al.ll2002ft . Given that 
there are more than 400 confirmed globular clusters within a 
radius of about 30 kpc dBarmbv et al.ll2000l) compared with 
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about 150 at all radii within the Milky Way, it is likely that 
there will be about 50 clusters orbiting M31 at radii useful 
for probing the halo. These clusters will be spread over a 
large area of the sky - spectroscopic confirmation of the na- 
ture of potential cluster candidates is essential, making this 
program very time-consuming from the ground. The accu- 
racy which the RVS will be able to achieve for integrated 
spectra has yet to be established. However, velocity errors 
of about 10—15 km s" 1 would be acceptable for constraining 
the halo mass of M31. RVS spectra can also be used to iden- 
tify other likely cluster candidates which can subsequently 
be followed-up from the ground. 

Bright (I ~ 16) AGB stars belonging to the halo and 
intermediate-age disc populations will provide additional 
tracers of the M31 potential. The ~ 10 4 young disc stars 
brighter than I = 16 and located throughout the M31 disc 
will permit a detailed comparison of th e stellar kine matics 
with the well-studied gas kinematics fe.g lBraunll99 lT) . They 
will also facilitate the detailed modelling of disc features such 
as the warp. 

2.16 Serendipitous discoveries 

Given the richness of the Gaia data set, it is to be expected 
that it will also throw up many surprises in the field of 
galactic dynamics. In addition to revolutionising our exist- 
ing picture of dynamical structures such as open clusters, 
moving groups, etc., many new features will undoubtedly 
be found. The scope for such discoveries is highlighted by 
a number of recent discoveries relating to the Milky Way 
disc, namely (i) the detection of a population of thick disc 
stars whose surprising kinematics suggest that they may 
have originated in the satellite galaxy whose merger with 
the original Milky Way d isc led to the formation of the thick 
disc iGilmore et all2002j) ; (ii) the identification of a ring-like 
struc ture at the outer edge of the Galactic disc lllbata et alJ 
120031) which spans roughly 100° in longitude and ±30° in lat- 
itude; (iii) the detection of an overdensity in the 2MASS all- 
sky M-giant catalog ue which may be th e remnant of an ac- 
creted dwarf galaxy iMartin et al.ll2004T) . The detailed rela- 
tionship between these observations is unclear. iHelmi et alJ 
J2003al) discuss a de bris origin for the r ing structure which 
has been shown by lYannv et alJ ll2003l) to be a kinemati- 
cally coherent structure. It has also been sugg ested that the 
overdensity identified bv lMart in et alJ j2004f) was the pro- 
genitor of the ring - howeve7 HM^m^mve^a!Tll2004l) suggest 
that the overdensity is instead related to the warp and flare 
in the external disc. However, what is clear is that the disc 
of the Milky Way has had a complicated formation history 
which Gaia will be uniquely able to probe. 

While these discoveries rely on observations which 
probe to magnitudes fainter than the cut-off for the RVS, it 
is important to note that the identification of the outer ring 
relied on the availability of large area photometric surveys: 
earlier, small-area surveys had simply overlooked its pres- 
ence due to its relatively low surface density. The true nature 
of this structure will only be confirmed once spectroscopy is 
available for more of its members in order to determine its 
kinematic relationship to the rest of the Milky Way disc. 
The discovery of the new thick disc population was based 
on radial velocities with accuracies of about 15kms _ , a 
level of accuracy which will be achieved by the RVS for all 
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Figure 10. Example spectra of common classes of outburst ob- 
jects for which Gaia will issue alerts to permit ground-based 
follow-up. See text for a detailed discussion. 



the stars it surveys. The great strength of the RVS is that 
any new structures which are identified will already have a 
wealth of kinematic data with which to determine their na- 
ture. If necessary, follow-up observations from the ground 
can be used to study the detailed properties of individual 
structures. 

During its five year lifetime, the Gaia satellite will ob- 
serve many transi ent objects, such as supe rnovae and mi- 
crolensing events ijBelokurov fc EvanslEooll . It is intended 
that many such events will be identified in real time on 
the ground, with alerts being issued (within 24 hours in 
the case of supernovae) to enable ground-based follow-up. 
RVS spectra will provide immediate classification of the 
brightest photometric alerts, those of greatest interest for 
follow-up. All photometric alerts with Cousin's 7c<14 mag 
will be bright enough for Gaia spectroscopy to provide dis- 
criminant classification on even the single, first epoch spec- 
trum. Fig. 1101 shows ground-based spectra, recorded with 
the Asiago 1.82m telescope operating in Gaia- like mode, of 
objects most frequently appearing among those undergoing 
outbursts. These spectra outline the major spectroscopic dif- 
ferences to be expected between different classes of outburst 
which will trigger Gaia photometric alerts. Brightening in 
pre-ZAMS objects (like T Tau and AB Aur in Fig. EH can 
lead to an increase of some magnitudes, accompanied by 
strong and wide Call emission (and weaker Paschen lines) 
from the circumstellar disc. The width of the emission lines 
in novae (Nova Cyg 2001-1 in Fig. HOB is very large, tracing 
the outward velocity of the ejected material (of the order of 
1 000 km sec -1 ). Dwarf novae (SS Cyg in Fig. 1101 display 
a smaller width connected with the lower Keplerian veloc- 
ities in the accretion disc, and the emission lines in symbi- 
otic novae (PU Vul in Fig. 1101 are even narrower tracing 
the small velocities of the wind from a late type giant sud- 
denly ionised by the burst of hard radiation from a white 
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dwarf (WD) companion undergoing a non-degenerate ther- 
monuclear outburst. The spectra of supernovae (type II SN 
2004et in Fig. E3 appear as a flat continuum, the width of 
the entire Gaia wavelength range corresponding to a typical 
ejection velocity of 10 000 km sec -1 . Finally, ejection of opti- 
cally thick material gives rise to spectacular P-Cyg profiles, 
like those displayed by V838 Mon (500 km sec -1 terminal 
velocity for the spectrum in Fig. 1101 later to became fa- 
mous for its huge light-echo. As the Figure clearly shows, 
RVS spectra will be able to distinguish between the differ- 
ent types of outburst thereby facilitating decision-making 
about the necessity and urgency of ground-based follow-up. 



3 BINARY STARS 

Another domain in which the Gaia mission will bring new 
and extensive data is the detection, classification and com- 
plete characterisation of binary systems. In this field, the 
RVS will be a key contributor by providing multi-epoch mea- 
surements of the radial velocity. It will impact on nearly all 
aspects of astrophysics where binarity is a basic quantity or 
provides the means to derive fundamental parameters. The 
statistical properties of binary stars (binary fraction, mass 
ratio and period distribution) will improve our understand- 
ing of the conditions of star formation in various locations in 
the Galaxy. Masses and radii, in addition to absolute lumi- 
nosities, are fundamental inputs to stellar physics (sec §0J). 
A knowledge of masses for a large variety of stellar types will 
improve the determination of the mass- luminosity function 
and directly impact on our knowledge of the Initial Mass 
Function (see § I3.2|l . 

The difficulty of obtaining unbiased samples from 
ground-based observations is underlined by the signifi- 
can t temporal separation of IDuauennov fc Mavorl (ll99lT) 
and lHalbwachs et all Jgffl ^ despite being papers II and III 
of the same series. In lDuauennov fc Mavorl l)l99lfl . 37 spec- 
troscopic orbits were derived and the duplicity of Solar-like 
stars in the Sola r neighbourhood was stu died using a sample 
of 164 systems. lHalbwachs et all l)2003h . which is based on 
two CORAVEL radial velocity surveys and Hipparcos dis- 
tances, uses an original sample of about 600 stars. The lim- 
itation of both papers comes from small number statistics. 
For comparison, Gaia will discover about 60 million binaries: 
more than 10 7 astrometric binaries, 10 6 eclipsing binaries, 
10 6 spectroscopic binaries, 10 7 resolved binaries within 250 
pc. Of these, some 10 000 will yield s tellar mass determ i- 
nations accurate to within 1 per cent jArenou et alj|2002>) . 
More specifically, the RVS will be the only instrument on 
board Gaia which will be able to discover virtually all non- 
eclipsing binaries with orb ital periods up to a mo nth for 
stars brighter than V = 14 dZwitter fc Munarill2004l) . These 
observations will provide large and unbiased samples of bi- 
naries for a very large range of periods and mass ratios and 
for stars at all evolutionary stages (even the shortest-lived). 
This will make it possible to study period and mass ratio 
distributions as a function of spectral type, population, for- 
mation site, etc. 

The specific cases of eclipsing binary stars fS 13.11 . cat- 
aclysmic variables f§ 13.31 . and symbiotic stars (S 13.41 are 
detailed below. 



3.1 Eclipsing Binary Stars 

In the case of eclipsing systems the radial velocity informa- 
tion will be combined with multi-band photometry to de- 
rive a complete orbital solution. Note that Gaia photometry 
alone can be used to discover much fainter binaries, but it 
will be limited to eclipsing systems and these are uncommon 
for periods longer than a few days. In addition, the RVS ve- 
locities will efficiently complem ent the astrometr ic data for 
the period range O.lyr to 5yr (ISoderhiermll2003ft . The ob- 
servation of eclipsing binaries by Gaia-RVS will be of prime 
importance to advances in stellar astrophysics, since stellar 
masses, in addition to radii and surface temperatures, can 
be derived with high accuracy. The distribution of param- 
eters such as mass ratios, eccentricity and orbital periods 
can cast light on the various processes which hav e been pro- 
posed for the formation of close binary systems feate et alJ 
2002). In addition, there are aspects of the binary evolu- 
tion models that are not well understood. We note that an 
apparent discrepancy between observations and models for 
binary masses near 0.8Mq has been found JPopperlll997al : 
IClausen et~al"ll200lf) . In these systems, the secondary star 
appears to be more evolved than the primary in both lu- 
minosity and radius. Modelling these stars is further com- 
plicated by the fact that many objects show evidence of 
spot activity. Good samples are needed before drawing any 
firm conclusions. The high resolution of the RVS spectra will 
allow the derivation of the statistical properties of spectro- 
scopic binaries with an accuracy that cannot be achieved 
from ground-based surveys, due to the required number of 
measurements. 

Gaia will der ive stellar masses to an accuracy better 
than 2 per cent llMarrese et al.l 120041 : IZwitter et all 120031 : 
iMunari et al.ll200lf) . stellar radii with 1 — 4 per cent accu- 
racy, and mass ratios with 1 per cent accuracy for more 
than 10 eclipsing binaries brighter than V=15; of those 
at least 25 per cent will also be double-lined spectroscopic 
binaries, providing an enormous data base. In addition, syn- 
thetic modelling of the spectra of spectroscopic binaries us- 
ing lines from both objects allows one to derive fundamental 
parameters such as temperature and surface gravity which 
can be directly compared with the results of the orbital so- 
lution. Those data can be directly compared with stellar 
models on the theoretical L-T e ff plane, thereby overcoming 
the difficulties (for example, bolometric corrections, colour- 
temperature transformations, distance, reddening) that al- 
ways hamper such a comparison on the observational V, (B- 
V) plane. If the stars of a binary pair have significantly dif- 
ferent masses, then the requirement that the models fit the 
data for a single age provides a strong constraint. An accu- 
racy of around 1 per cent in mass and radius is then required 
to retrieve accurate informa tion about opaci ties, convection 
prescriptions and rotation (lAnderser]|2002D . The need for 
additional mixing in stars has been confirmed by several 
authors using binary s tar data. We note, am ong others, 
th e pioneering wor k by lAndersen et alJ lll988l). the review 
bv I Andersen! I2OO2T1 , and the recent work bv lSiviero et alJ 
i2004f) who find excellent agreement between the Padova 
isochrones and the obse r ved d ata on both components of 
V432 Aur. lYoung et alJ (1200111 present new evidence that 
the overshoot efficiency might depend on the stellar mass. 
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All these results need to be tested on the very much larger 
data base that Gaia will provide. 

Knowledge of the temperatures and radii of stars accu- 
rately determines their luminosities. Combining these with 
measured mul ti-colour apparent m agnitudes and allowing 
for reddening llPrsa fc Zwitterl l200Bl) yields an accurate dis- 
tance to the binary JWilson fc Wvithel2003l) . Note that this 
method of distance determination is complementary to as- 
trometric measurements and is limited only by the limiting 
magnitude of the instrument. The meth od has recently been 
used to derive a distance to the Pleiades ijMunari et all2f)o3 l 
and to study the dista nces of objects as far a s the Large 
Magellanic Cloud (see iFitzpatrick et all 120031 . and refer- 
ences therein) . Results of the binary method agree well with 
astrometric distan ces for a dozen eclipsing binaries discov- 
ered by Hipparcos jZwitter fc Munarll2004l) . 

3.2 Mass Function Determination 

The stellar initial mass function (IMF) is of fundamental 
importance in many fields. On the one hand, the stellar 
mass distribution determines the evolution, surface bright- 
ness, chemical enrichment, and baryonic content of galaxies. 
On the other hand, knowledge of the slope of the IMF and its 
possible variations can cast light on the physics of cloud frag- 
mentation and the star formation processes through which 
different mass ranges of star are assembled. 

Determining the IMF of a stellar population of mixed 
ages is a cumbersome affair. In fact, the IMF is usually ob- 
tained from the stellar luminosity function (SLF), i.e. the 
number of stars in a survey volume per magnitude interval, 
through knowledge of the mass-luminosity relation (MLR) 
in all mass ranges and of the star formation rate. In the 
vast majority of cases, the stellar masses cannot be directly 
derived and the mass has to be deduced indirectly from the 
luminosity and the evolutionary state of a star. Two basic 
approaches have been tried in the literature to derive the 
SLF: the first one makes use of a local volume-limited cat- 
alogue of stars with well-measured distances. The second 
method takes larger samples of stars from deep photome- 
try. This latter is affected by various spurious effects, such 
as the Malmquist bias, imprecise determi nation of the c om- 
pleteness a nd unknown bi nary corrections <Kroupal200llL In 
particular. iKroupal lll995l) has shown that the significant dif- 
ference between the volume-limited and magnitude-limited 
SLF for the local disc population at magnitudes fainter than 
My = 11.5 is mainly due to the presence of unresolved bi- 
nary stars. 

The determination of the SLF from a volume-limited 
star catalogue gives more reliable results. However, to date, 
complete samples of trigonometric parallaxes are known only 
for stars brighter than My ~ 9.5 at distances d < 20 pc, 
while for the f aint M dwarfs the es timated completeness dis- 
tance is 5 pc (iLeinert et al.lll99^ . This means that only a 
limited sample of objects is covered. A major caveat of any 
photometric LF is that the determination of the distance re- 
lies on the photometric determination from a CMD. In prac- 
tice, the determination of the IMF from the SLF requires 
knowledge of the chemical composition of each star, since 
the absolute magnitude and colour depend on both metal- 
licity and age. Assuming Solar metallicity for a metal poor 
thick disc star would lead to an underestimate of the abso- 



lute magnitude, an overestimate of the distance and thus an 
underestimate of the number density. 

The RVS will contribute to the determination of the 
IMF through the accurate estimation of the binary fraction 
and of stellar parameters. This, together with a parallax de- 
termination, will precisely locate a star on the CMD and will 
allow a better determination of the stellar MLR which can 
subsequently be used to convert stellar luminosity functions 
into stellar mass functions. Our present-day knowledge of 
the MLR for masses lower than 1.5Mq is based on data on 
visual binaries, since only a limited amount of high quality 
data on double-lined eclipsing binarie s and resolved spect ro- 
scopic binaries is currently available llMalkov et al]|l997l) . 

By combining the proper motion, radial velocity infor- 
mation and chemical abundances of single stars, information 
can be derived about the dependence of the IMF on metal 
content. In fact, while at the upper mass end (> lOM©) the 
IMF seems to be virtually independent of the metal con- 
tent, a dependence ha s been suggested at the low mass end 
jRevle fc Robinll200ll) where a flatter slope has been found 
for thick disc stars. However, while on the one hand the 
IMF of the metal poor st ars in the Galactic spheroid does 
not s how this behaviour iGould et allll998fl . on the other 
hand lZoccali et all feOOOaf) derive an IMF for the metal rich 
Galactic bulge which is consistent with that of the metal 
poor globular clusters. 

The existence, if any, of extremely metal poor stars will 
cast light on the primordial IMF. Any observed variation of 
the IMF between different environments would suggest that 
cloud fragmentation based on the Jeans formulation (i.e. on 
gravity) is not the only mechanism at work, but rather that 
different processes (such as turbulence) m ight be important 
jKlessenll200llPadoan fc Nordhmdlll999Tl . 

3.3 Cataclysmic Variables 

Cataclysmic Variables (CVs) are a diverse class of short- 
period semi-detached binaries consisting of an accreting 
white dwarf (WD) primary and (typically) a low-mass, 
main-sequence secondary star. They are valuable manifes- 
tations of late-stage binary evolution, which also provide 
a window into the fundamental physical processes associ- 
ated with accretion and nova explosions (including, possi- 
bly, Type la supernovae). The RVS spectra will be useful 
for identifying CVs and will contain substantial astrophysi- 
cal information on line strengths, shapes and velocity vari- 
ations. CV spectra show strong and distinctive lines in the 
RVS band: broad (~ 3000 km s -1 , 8nm) double-peaked Ca 
lines, originating in the cooler outer edges of the disc; nar- 
row emission (from the heated face of the secondary star) 
or, in longer period systems where the secondary is larger 
and brighter in this band, photospheric absorption features 
instead of emission; broad Paschen lines, one of which coin- 
cides with the Ca line at ~ 866 nm. The spectra vary sub- 
stantially from system to system, depending on the orbital 
period, accretion rate and disc state. Uniquely among as- 
trophysical objects, some of these systems show extremely 
strong NI lines. In magnetic CVs, the Ca triplet can be seen 
strongly in emission in the one or two systems studied but 
very few data are available. 

For the current RVS specifications a practical limiting 
magnitude of V ~ 16 can be expected for determining CV 



© 0000 RAS, MNRAS 000, 000-000 



Expected science yield from the Gaia Radial Velocity Spectrometer 23 



system para meters. There are cu rrently 140 CVs brighter 
than V=16 iDownes et al.ll200lT) . Extrapola ting from the 
initia l findings of the Sloan Digital Sky Survey l|Szkodv et alJ 
120021) . the incompleteness level is about 30 per cent on the 
basis of selection by colour only (which has strong selec- 
tion effects). Gaia spectroscopy will be very sensitive in the 
search for signatures of accretion (line emission). This should 
lead to the discovery of a large number of intrinsically faint 
systems (possibly up to 1000). All of these CVs will have 
excellent parallaxes and hence luminosities. The ratios of 
secondary to primary masses provide an important test of 
population models of interacting binaries and can be calcu- 
lated from the orbital period and velocity of the secondary: 
for a significant fraction of CVs in the RVS sample, the bi- 
nary orbit and inclination will be obtained from reflection 
modelling, eclipses, or directly from the Gaia astrometry. 

Gaia will provide a minimally biased sample of CVs 
for population and evolution studies. It should be noted 
that spectroscopy from the RVS is particularly powerful: 
many CVs hardly vary (the nova-likes) and many systems 
which would be considered detached based on photometry 
are evidently accreting only once spectroscopy is obtained. 
Periods are much easier to determine with radial veloc- 
ity measurements than by photometry. The spectroscopy 
will also pick up new classes of unexpected objects such 
as short period systems with K- rather than M-type secon- 
daries. There is also the likelihood of finding longer period 
double-degenerate and other 'graveyard' CVs, as evolution 
models predict large numbers of these. In current surveys 
they would be indistinguishable from WDs on the basis of 
their colours. In addition, the RVS survey will determine 
the fraction with Ca emission compared to absorption and 
the relationship of these with the different classes of CV 
(magnetic/non-magnetic), luminosity and secondary spec- 
tral type. Ultimately, information will be gained on the cause 
of the differing line strengths. 

The Ca triplet lines are more clearly double-peaked 
than other strong lines in the optical spectrum because they 
originate in the cooler outer regions of the disc. For a large 
number of disc CVs, Doppler tomography will be possible. 
This will produce maps of disc velocities including the ef- 
fects of tidal distortions and the infalling stream, and the 
run of Ca emission within the disc. It will also identify other 
aspects of disc structure, for example the spiral waves tha t 
have been seen in some Dwarf Novae dSteeehs et alJll997ft . 
which have implications for our understanding of accretion 
discs in general. For those CVs for which data are avail- 
able, many (perhaps most) seem to have narrow Ca triplet 
emission components. The strength of these lines and their 
ratios are inputs for atmospheric heating models of the sec- 
ondary star. The Ca triplet is particularly good for both of 
the above applications: Ca velocity maps are of intrinsically 
higher resolution than the more commonly used Ha maps 
and, because Ca is less saturated than Ha, the associated 
atmo spheric heating models are simplified iMarsh fc Duck! 
11996ft . 

3.4 Symbiotic stars: the Supernova connection 

Type la Supernovae (SN la), the only supernovae to ex- 
plode in elliptical galaxies, are widely used as cosmologi- 
cal distance indicators. They provide a clear indication of 



cosmic acceleration and of a non-zero value for E instein's 
cosmological constant (e.g. iPerlmutter et af 1 11999ft . While 
current search programs have been successful in discover- 
ing new SN la at high redshift (e.g. the High-Z Supernova 
Search Team (HZSST) and Supernova Cosmology Project 
(SCP) surveys), the true nature of the SN la progenitors 
and their explosion mechanisms is still a matter of inves- 
tigation and debate. With look-back times of roughly half 
the current age of the Universe, one has to make sure that 
possible evolution of the progenitors is not mimicking a cos- 
mological effect. It is generally believed that the explosion 
of a SN la is associated with a WD that grows in mass 
to the Chandrasekhar limit via accretion. Most of the de- 
bate about the SN la progenitors involves the nature of the 
companion feeding mass to the WD. In the single degener- 
ate scenario it is a main-sequence or giant star, while in the 
double degenerate scenario the companion is itself a WD and 
the supernova explosion results from the merger of the two 
WDs in a binary system. 

A possibility for the single degenerate scenario is that 
the donor star is a late type g i ant, giving rise to a sym- 
biotic star (SyS) <KenvoiJil98al . iMunarl (Il994l) has listed 
the requirements that SyS have to possess to be viable pro- 
genitors of SN la. These include: (a) their WDs must be 
massive enough to be able to grow to the Chandrasekhar 
limit by accretion of only a fraction of the mass reservoir 
transferable from their companions; (b) their number in the 
Galaxy must be appropriate to account for the SN la rate 
given the evolutionary time scale of the donor giant; (c) the 
SyS must be appropriately frequent among the bulge/thick- 
disc/halo populations of the Galaxy, which correspond to 
the dominant stellar populations in elliptical galaxies. 

Gaia is ideally suited to address all these issues. First, 
Gaia will detect most late type giants in the Galaxy and 
it will be able to discover most of those harboured in 
SyS systems, thus providing the fractional rates within 
the various Galactic populations. SyS systems display a 
distinctive appearance in the RVS wavelength range (e.g. 
Schulte-LadbeclJll988t IZhu et alJ Il999l: iMunari fc Zwitterl 
20021: iMarrese et alJl20o3T that, coupled with their charac- 
teristic photometric variability and spectral energy distri- 
bution, will allow the RVS to identify them easily even in 
low S/N spectra. In fact, the cool giants in Galactic SyS are 
mainly of spectral types K and M, thus emitting strongly 
in the Call region, away from the veiling effect of the blue 
continuum emitted by the circumstellar ionised gas. This 
gas also shows up in the RVS spectral range by display- 
ing strong Paschen, Call and Hel emission lines (sometimes 
also useful to trace the orbital motion of the ionising WD 
companion). Monitoring their radial velocities over the 5 yr 
mission life-time (3 times longer than the typical ~ 600 day 
orbital periods), the RVS will provide orbital periods for a 
statistically significant number of objects. An indication of 
the mass ratio between the cool giant and the accreting WD 
can be obtained by comparing the radial velocity curve of 
the giant with that of the emission lines of the highest ion- 
isation states which presumably arise in the regions closest 
to the WD. In fact, none of the strong lines of very high 
ionisation levels in the spectra of SyS fall within the RVS 
wavelength range. However, knowing their orbital periods 
and phases from the combination of Gaia photometry and 
radial velocities, a small number of ground-based spectra 
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at appropriate wavelengths and phases will suffice to deter- 
mine the mass ratio for the SyS discovered by Gaia. These 
will then constrain the mass distribution function of WDs 
in SyS. The Gaia medium-band photometric bands in the 
red part of the spectrum have been designed to be sensitive 
to molecular bands. The intensity of these bands is metallic- 
ity dependent, which yields an additional constraint on the 
metallicity of the SyS. Finally, the location of the SyS in 
the Galactic phase-space, coupled with the metallicity and 
chemical analysis, will firmly address their relationship to 
the various Galactic stellar populations and thereby con- 
stra in their rate of occurrenc e in elliptical galaxies. 

iMunari fc Renziml l|l992f ) estimated the number of SyS 
in the whole Galaxy to be about 3 x 10° , which implies that if 
only 2 — 4 per cent of them end their evolution in a supernova 
explosion it would acc ount for the observed SN la rate. Their 
Galactic kinematics jMunari fc Busonl Il994h and infrared 
photometric properties ijWhitelock fc Munaril Il992h argue 
that the majority of SyS are associated with the bulge/thick- 
disc component of the Galaxy. Evidence from energy radi- 
ated in the UV l|Munari fc Busonl ITifflih and from orbital 
motions suggest that a sizeable fraction of known SyS har- 
bour massive WDs, resembling stable H-burning conditions 
and only minor mass loss. Thus, all the ingredients seem to 
be in place but it will be Gaia that will provide the statis- 
tically sufficient data set to address definitively the issue of 
whether or not SyS are viable SN la progenitors. 



4 TOWARDS AN UNDERSTANDING OF 
STELLAR EVOLUTION 

4.1 Introduction: RVS impact on stellar 
astrophysics 

In spite of considerable recent efforts in the area of stel- 
lar evolution, many open problems still remain. The micro- 
physics of the equation of state and opacities has been 
tested in the range of parameters corresponding to the ra- 
diative zone of the Sun using asteroseismology. However, our 
knowledge of macro-physics processes such as rotation, con- 
vection and turbulence over the entire stellar mass range 
remains very poor. In order to improve significantly our 
knowledge of stellar interiors, it is necessary to couple the 
determination of global parameters for a statistically sig- 
nificant sample of stars with information from seismology 
ijLebreton fc Baglinll2002l : lLebretorll200l l200oh . Seismology 
needs a precise knowledge of the global parameters of the 
stars: luminosity, effective temperature, radii or masses in 
the case of binaries, surface chemical composition. In fact, 
seismology can be used to derive the fundam ental stellar pa- 
rameters from stellar oscillation frequen cies iPopperlll997tJ : 
iPetersen fc Christensen- Dalsgaardlll996ft . However, if global 
parameters are known, stellar oscillation frequencies can be 
predicted and compared with the observations, giving infor- 
mation about the physics of the stellar interiors. 

Gaia astrometry, photometry and spectroscopy will 
build a complete and homogeneous sample of accurate global 
parameters for a large range of stellar masses. The spec- 
tra from the RVS in particular can be used to derive the 
effective temperature T e ff, surface gravity logg, metallic- 
ity [Fe/H], rotational velocity V ro t, chromospheric activity 



and interstellar extinction (see Tableland Paper I). When 
combined with parallax determ inations, the lumin osity, ra- 
dius and mass can be obtained jBailer-Jonesll2005D . The ex- 
pected accuracies of the various parameter determinations 
are discussed in Paper I (see also Table 0. Finally, the ab- 
solute magnitudes of stars having distances measured with 
an accuracy of 1 per cent will be known to within 0.03 mag. 
The accurate placement of stars in HR diagrams will allow 
the construction of very precise stellar tracks and give hints 
about the physical processes taking place in stellar interiors. 

It is well known that the location of a star in the HR di- 
agram does not allow a unique determination of its age, since 
several combinations of [Fe/H], [a/Fe] and age are possible. 
To derive accurate ages requires very precise estimates of ef- 
fective temperatures, reddening and chemical composition. 
For old main sequence halo stars the expected RVS accu- 
racy on the stellar parameters (see Table 0, combined with 
information from the photometric observations, will lead to 
uncertainties in individual age determinations of about 20 
per cent. The largest contribution to this error is from the 
uncertainty in [Fe/H]. However, since the age determination 
will be made by comparing the location of the stars in the 
HR diagram with theoretical models, further uncertainties 
will come from our poor knowledge of stellar models (i.e. 
mixing, diffusion, nuclear rates, etc.). The stellar parameter 
determinations from Gaia will be therefore first be used as 
input to improve our knowledge of stellar models. 

In the following sections, we underline some of the ar- 
eas where the RVS will particularly contribute to our under- 
standing of stellar physics. 

4.2 Star formation and pre-main-sequence 
evolution 

The process of star formation in different environments is far 
from being understood. In particular, it is difficult to recon- 
cile the prominent influence of the local environment (tur- 
bulence, compression, initial trigger) on small scales with 
the universalit y of the Schmidt and Kennicut law on Galac- 
tic scales (see lElmegreenll200l for a recent review) which 
suggests that Galactic-scale gravity is involved in the first 
stages of star formation. The Gaia study of Galactic star 
forming regions will shed light on the details and modal- 
ity of star formation both in space (does the star formation 
take place in giant molecular clouds or in small cloudlets?), 
and in time (from the age distribution of stars) deriving 
absolute luminosities, effective temperatures, multiplicities, 
kinematics and, where possible, masses of the star forming 
complexes. 

Many aspects need to be clarified concerning the pre- 
main-sequence (PMS) evolution of stars. In the domain of 
small masses (< I.5M0) the main uncertainties are related 
to the treatment of over-adiabatic convection, the stars be- 
ing fully convective at the beginning and developing a small 
radiative core during the H-burning phase. The PMS tracks 
run almost parallel to the main sequence, starting from the 
deuterium burning phase which represents the starting point 
of the evolution in the visible. Unless binaries are present, 
mass e stimates rely on the assumptions of the convective 
model fo'Antonalll999l) . The prototypical PMS star of this 
mass range is T Tauri. During those initial stages, param- 
eters such as stellar rotation, angular momentum evolution 
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and magnetic field strength play very important roles which 
need to be quantified by observed data. 

Massive stars are believed to be for med via co ntinu- 
ous mass accretion during the PMS phase (Palla 1998). The 
star follows the birth-line accreting mass and, when the 
accretion phase stops, the object leaves the birth-line and 
moves towards the main sequence. This model is in agree- 
ment with several observational constraints. However, not 
everything is well explained. The location of the birth-line 
is strongly dependent on the accretion rate which is essen- 
tially unknown. The accretion rate greatly influences all the 
internal properties of the stars as well as the PMS lifetime 
iMaeder fc Behrendll2002T) . Additionally, the accretion rate 
derived for low-mass stars fails to describe massive objects, 
since the formation time would be too long compa red to 
the main- sequence lifetime (iNorberg fc Maederllioool) : stars 
would leave the main sequence before being fully formed. 

An additional source of uncertainty for massive PMS 
stars is related to their high rotation velocities and to the 
treatment of convection (see the following Sections for a de- 
tailed discussion). Accurate distance determinations for very 
young associations and clusters, together with spectroscopic 
determination of the effective temperatures, rotational ve- 
locities and infrared spectra of PMS stars will lead to more 
accurate PMS tracks. 

The angular momentum budget of stars and the way 
in which stars lose it has not yet been fully understood. 
PMS stars are known to be fast rotators. In young clus- 
ters there seems to be a relation between the ages and ro- 
tational velocities of the stars, in th e sense that stars ap- 
pear to be losing angular m omentum jTerndrup et alj|2000t 
IStassun fc TerndrudEioll) during their lifetime. The trans- 
fer of angular momentum through the star and the role of 
the proto-stellar disc predicted by the models need to be 
compared with observations. 

A large number of (proto)-stars in the PMS phase are 
likely to be detected by Gaia, with a significant sample of all 
ages ranging from the earliest T-Tauri phase (few Myr) to 
the onset of central H-burning (several tens of Myr). Gaia 
will observe about 120 young open clusters like Praesepe 
within 1 kpc of the Sun with an accuracy comparable to 
that reached by Hipparcos in the Hyades where low mass 
PMS stars are visible. In the RVS wavelength range, the 
spectra of T-Tauri stars show all three Call lines in very 
strong emission. This will allow for the identification of T- 
Tauri stars even at very low S/N ratios. 

4.3 Rotation 

Despite the fact that many observations have demonstrated 
that r otation is a necessary ingredient in models of ma ssive 
stars dMaeder fc MevnetfeOOObl ISoderblom et al.ll200lf) . the 
effect of stellar rotation on the evolution of stars has only 
recently been included in stellar models. Among the critical 
observations, we recall the fact that rapidly rotating massive 
O stars have peculiar He abundances. B- and A-type stars 
in the Magellanic Clouds are found to have large relative ex- 
cesses of He/H with respect to the prediction of current mod- 
els without rotation. The most likely explanation of these pe- 
culiar abundances is rotation-induced mixing. Rotation can 
influence the location of a star in the HR diagram, its lumi- 
nosity and its lifetime. The evolution of massive stars is the 



result of an interplay between rapid rotation and mass loss. 
In fact, massive objects are rapidly rotating when they form 
and as a result of the reduced internal pressure they are sub- 
luminous. During the main-sequence lifetime the stars decel- 
erate owing to angular momentum loss through their winds 
and b ecome more luminous mor e rapidly than non-rotating 
stars jMaeder fc MevnetlEoOOrl . Rotation increases the H- 
burning time for stars more massive than 9 Mq by about 
25—30 per cent, while the e ffect on the He-burning lifet ime is 
smaller (about 10 per cent: lMaeder fc Mevnetll2000bl) . The 
effect on age determinations is far from being negligible, 
amounting to about 25 per cent. In that sense, rotation can 
mimic the effects of core overshoot. To date, only 20 00 
stars have measured vs'mi JClebocki fc Stawikowskll2000D . 
A detailed analysis has revealed that while early type stars, 
from O-type to early F, are fast rotators (50 — 400 kms -1 ), 
late type objects (from late F to M ) possess low e r rota- 
tional velocities ( usini < 50 kms -1 : iMunari et alj . l2001al . 
ISoderblom et all 1200 if) . PMS stars are known to be fast 
rotators reaching 30 per cent of the breakup velocity - as 
they age, mass loss causes them to lose angular momentum. 
In fact, in young star clusters, a relation has been found be- 
tween the age of the cl uster and the distribution of rotational 
velocities of the stars jTerndrup et all200ll20o3) . Gaia will 
study 120 young clusters closer than 1 kpc, allowing the de- 
termination of the rotational velocity. The interaction of the 
convective envelope of low mass stars with their differential 
rotation sustains a magnetic field which, by trapping the 
ionised stellar winds, induces a loss of angular momentum 
from the stars. 

The RVS will greatly improve our knowledge by directly 
measuring rotation via line broadening and spot transit for a 
large sample of stars: the combination of accurate luminos- 
ity, effective temperature and rotational velocity determina- 
tions prov ided by the RV S sp ectra will cast light on their 
evolution. ICombod J2003t) and lGomboc fc Katzl J2005I) dis- 
cuss the accuracy of the determination of usini from the 
cross-correlation of the RVS spectra of single stars. Preci- 
sions of wsini ~ 5 km/s should be obtained at the end of 
the mission for late type stars down to V=15. For B5 main- 
sequence stars, we expect an accuracy of v sin i ~ 10-20 km/s 
up to V~ 10-11. The periodic transit of spots on the stellar 
surface can also trace the rotation. The RVS spectra can 
trace several spots at once as they cross the projected sur- 
face. The associated emission lines will split according to the 
velocity, allowing accurate rotation measu rements even with 
only a few tens of spectra jMunarill2003rl . 

4.4 Mixing processes 

One of the main points of uncertainty in stellar evolution 
theory concerns the extent of the convectively unstable re- 
gions and associated mixing. Thermal convection arises in 
stars when radiation is not sufficient to carry the heat flux 
originating from the deep interior. The convective insta- 
bility occurs wherever the local temperature gradient is 
steeper than the adiabatic gradient, a condition called the 
Schwarzschild criterion. 

However, an extra mixing beyond the classical con- 
vective regio ns (overshoot) seems to be at work (see 
IChiosil Il999l . for a review). The theory of non-local 
convection has made significant progress over the past 
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decad e JXiond ll99Ct ICanutd fl999l: ICanuto fc Dubovikoyl 
19981: ICanuto fc MazziteUi Il992j: iGrossmanl Il996t 
Ventura et all Il998l: ICanutd Il999l. hood : iBrummell et alJ 
20021: lYoung fc Arnetdl2002l) . However, the lack of satisfac- 
tory results means that the majority of stellar models are 
still calculated using the local approa ch of m i xing l ength 
theory (from t he pi oneering work by iMaeder j Jl975l) and 
iBressan et alJ il98lf) to the recent work of lGiraTdi et alJ 
l2002T) T. where the mean free path of the convective ele- 
ments is proportional to the scale height of the pressure H p 
through a free parameter derived from comparisons with 
the data. Overshoot can greatly influence the evolution of 
massive and intermediate-mass stars, changing the lifetimes 
in the H- and He-burning phases, as well as the luminosity 
of the stars. The effect on the age determination of a stellar 
population is far from being negligible and is about 25 per 
cent for an A star of age 2 Gyr. Element diffusion can 
also change the evolution of a star, which can introduce 
errors as large as 100 Myr in the age of a 1.7M© star. As 
far as low mass stars are concerned, in the mass range 
where the core switches from the radiative to the convective 
regime the determination of the overshoot paramete r 
leads to unsatisfactory results dWoo fc Demarauel l200ll) . 
This corresponds to the age range of the oldest open 
clusters (6 Gyr), introducing large uncertainties on the 
age determination itself. By combining the expected Gaia 
determinations of the global parameters of stars with the 
results of asteroseismology measurements (see § 14. It . it has 
been estimated that when L and T e fi of stars are known 
with a precision of about 2 per cent, then the overshoot 
parameter can be derived with a precisio n of 0.03 times the 
pressure scale height H p l|Lebretonll200ft) . 

A discussion of the lines of o bservational e vidence for 
the presence of overshoot is given bv lChiosil<ll998ft . The anal- 
ysis of the HR diagram of Hipparcos field stars shows that 
the region of the so-called Hertzsprung gap is very sensitive 
to overs hoot on the main sequence for stars with masses of 
1.6 Mq flSchroeded 1 1998F) . Using th e data to infer the his - 
tory of star formation in the disc, iBertelli fc Nasil (|200ll) 
find that the compatibility between the number of main- 
sequence and post- main- sequence stars requires some extra 
mixing during the core H-burning phase in stars in the range 
1.5-2 Mq. Information about the size of the convective core 
can be derived from mHe the maximum mass undergoing a 
core He-flash. mHe can be derived from the luminosity func- 
tion of the red clump stars: the mass distribution inside the 
clump presents a peak at m He where the He-burning time 
in the core has a maximum JOirardi fc Sa larisl l200ll f. The 
ratio 12 C/ 13 C in AGB stars can be derived from RVS spec- 
troscopy and gives hints about the third dredge-up, which 
in turn depen ds on the mas s loss rate and on the treatment 
of convection IIMarigdl200Cl) . 

Classical Cepheids provide quantitative constraints on 
the efficiency of mixing phenomena (rotational mixing, 
overshoot, semi-convection ) among inter mediate-mass stars 
jMaeder fc Mevned EoQObl: ICassisil liooi) . The discrepancy 
between evolutionary masses and pulsation masses for Clas- 
sical Cepheids dates back to the 1970s. Recent investigations 
suggest that this discr epancy is at the level of 10 per cent 
for Galactic Cepheids (|Bono et al.ll20oih but of the order of 
15 — 20 per cent for Magellanic Cepheids feeaulieu et alJ 
l200ll) . We still lack clear physical arguments to explain 



whether this discrepancy is either a real feature or a con- 
sequence of errors in the physical assumptions adopted to 
const ruct evolutionary and pulsation models l)Bono et alJ 
120021) . Suggestions have been advanced in the literature that 
the Cephei d mass discrepancy can be alleviated by including 
overshoot ijChiosi et al.lll992lh . 

The careful calibration of the HR diagram as a function 
of age and metallicity by Gaia will give very important hints 
about stellar interiors and, when considered in conjunction 
with the results of seismology observations, will improve our 
capability to construct non-local convective models. 

4.5 Mass loss 

Mass loss from stars is important for the evolution of the 
interstellar medium whose composition depends on the na- 
ture of the supplied material. Stars are known to lose mass 
at various rates during different stages of their evolution. 
The RVS observations will be well adapted for the study of 
stellar winds and mass loss, since Call and H lines display 
characteristic P-Cyg profile s in the RVS sp ectral range even 
for modest mass loss rates jMunari|l2003bT) . 

In particul ar, mass loss is a domina nt effect in the upper 
main sequence ijChiosi fc Maederlll98fi ) . The mass loss rates 
currently used are based on o bservations dde Jaeer et alJ 
amcrs fc Cassinellilll996Tl . However, the relation be- 
tween rotation and mass loss is not yet properly quantified. 
While lVardval lll985l) finds a significant increase in mass loss 
with rot ation for OB stars, only a modest increase is de- 
rived by iNieuwenhuiizen fc de Jagerl lll988h . The same un- 
certainty is present on the theoretical side. Whether or not 
mass loss influences the angular momentum seems to depend 
on whether it is proceedin g via equatorial or polar winds (see 
iMaeder fc MevnedEoOOal. for a review). 

In the domain of low mass stars, it is not yet 
clear whether the pulsation instability might drive the 
efficiency of mass-loss close to the RGB tip. There is 
mounting empirical evidence that a substantial fraction 
of br ight RGB stars are variables fedmonds fc Gillilandl 
Il996h . This result was strongly supported by detailed anal- 
ysis of time-series data collected by microlens ing experi- 
ments jKiss fc Beddindl2003t IWood et al. II2003T) as well as 
by Near-Infrared surveys dlta et al.ll2002l) . This is a new 
and very interesting result because it has been suggested 
that these objects are excited stocha stically by convec- 
tion llChristensen-Palsgaard et al .1200 ll) . The unique oppor- 
tunity provided by Gaia to measure distances and to col- 
lect homogeneous multi-band photometric and spectroscopic 
data for a very large sample of halo and disc RGB stars is 
crucial to provide a comprehensive analysis of these objects. 
In particular, the RVS will supply firm constraints on the 
role played by chemical composition and binary companions 
in the pulsational de-stabilisation of K-type stars. 

4.6 Variable stars 

Variable stars are typically used as primary distance indica- 
tors or as stellar tracers in the Galaxy as well as in external 
galaxies. However, they also play a fundamental role in stel- 
lar astrophysics, because they can be easily identified and 
their radial velocity curves are key observables to estimate 
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fundamental parameters such as stellar masses, radii, and 
effective temperatures. Therefore, it is not surprising that 
variable stars are also adopted to constrain basic physics 
problems such as an upper limit to the neutrino magne tic 
moment iRaffehll990l:ICasterlani fc degl'Innocentll993l or 
the existence of dark matter particles " iJPeln^oni^tHiL^IQfl 
lRaffeltl2OO0l) . Moreover, they also provide the unique oppor- 
tunity to compare predictions based on hydrodynamical pul- 
sation models, evolutionary models, and stellar atmosphere 
models. Gaia will play a paramount role in the field of vari- 
able stars and th is topic has already been d i scussed in se v- 
eral papers (e.g. iTammann fc ReindllEx)! iBonol l2003allcl 
and references therein). Alongside trigonometric parallax 
determinations, Gaia will provide photometric and spectro- 
scopic time series data for a large number of variable stars 
throughout the Galactic spheroid. This will create a unique 
opportunity to address, on a quantitative basis, several long- 
standing problems. In the following, we briefly mention some 
of the key issues. 

4-6.1 Intermediate-mass Variable Stars 

Recent theoretical feono et alJ Il999alfbh and empiri- 
cal iTammann et al.ll2003l : iKervella et alJl2004) results in- 
dicate that the Cepheid Period-Luminosity (PL) relation 
is not universal. In fact, both the zero-point and the 
slope of the optical PL relation seem to depend on the 
metal content. Moreover, detailed investiga tions of Galac- 
tic and Magellanic fundamen tal Cepheids dSandage et alJ 
120041: iKanbur fe Njgow[J [200^ support theoretical predic- 
tions IIBTn^et^dTll999al lbll concerning the nonlinearity of 
optical PL relations when moving from short to long- 
periods. This has not been conclusively established, how- 
ever, since differ ent theoretical predictio ns based on linear 
Cepheid models fearaffe fc Alibertll200ll) and e mpirical ap- 
proaches based on the Infrared-Flux method (iLuck et alJ 
120031: IStorm et alJl200 j) suggest a mild dependence on the 
metal content. More recently, theoretical predictions based 
on nonlinear, convective models and multi-band observa- 
tions suggest that the PL relation of first-overtone Galac- 
tic and Ma gellanic Cepheids marginally depends on metal 
abundance iBono et alJl2002T) . 

Gaia will supply a new impetus to studies of the 
Cepheid distance scale because it will obtain an almost com- 
plete census of Galactic Cepheids. Classical Cepheids in the 
period range from 6 to 16 days present a well-defined bump 
in both their luminosity and radial velocity curves. This fea- 
ture, for periods shorter than about 9 days, is located along 
the decreasing branch, while for longer periods it is located 
along the rising branch. Fig. llll presents predicted radial ve- 
locity and light curves in different photometric bands for one 
such Bump Cepheid in the Large Magellanic Cloud. As the 
Figure illustrates, the velocity changes will be detectable by 
the RVS at a statistically significant level down to a mag- 
nitude limit of about V = 14, where the velocity errors are 
below 6kms _1 for single epoch observations (Paper I). Ac- 
curate chemical compositions (intrinsic error smaller than 
0.2 dex) will als o be obtained down to a limiting mag nitude 
of V « 12 - 13 llMunarill2003al: iThevenin et alJl2005f> . 

We note that the amplitude of radial velocity variations 
exhibited by Classical Cepheids ranges from 15 to 90kms _1 
as one moves from short period to long period objects (see 



e.g. lUdalski et al.ll999l : IPetroni et al.l2003t) . By assuming an 
absolute magnitude of My ~ — 2 for short-period Cepheids, 
the RVS will thus provide accurate measurements of chemi- 
cal composition, reddening, and pulsation properties for all 
Cepheids located in the outer disc out to distances of about 
6—10 kpc (depending on reddening) and over the entire 
Galactic disc for the brightest ones. 

The current observational data set of metal- 
licity measurements for Galactic and Magellanic 
Cepheids is quite limited. Accurate spectroscopic mea- 
surements of Galactic Cep h eids have been recently 
pro vided by | Frv fc Carney! (Il997l . 23 objects) and 
by lAndrievskv et al.l i2004 131 objects). But only for 
three dozen of these Cepheids are ac curate absolute dis - 
tance determinations also available dStorm et alJ l2004t) . 
Accurate spectroscopic measurements of M agellanic 
Cepheids are available for 40 objects jLuck et al.| |l9 9il; 
Romanicllo. Prima sTMottini. Groenewegen. Bono, fc Francoisl 



1200511 . but accurate distan ce determina t ions are available 
only for a few of them dGieren et al.l |2000F) . Moreover, 
optical PL relations based on nonlinear convective mod- 
els predict that, on moving from metal-poor (Z=0.004, 
the metal abundance typical of Small Magellanic Cloud 
Cepheids) to metal-rich (Z=0.02, the metal abundance 
typical of Galactic Cepheids), there is a difference in 
magnitude that at logP = 1 is 0.4 mag in the V — band and 
0.3 mag in the /—band jFiorentino et al.ll2002l) . Empirical 
estimates of this effect still present a broad range of values. 
Using the RVS data, it will be possible to perform a robust 
calibration of the optical PL and Period-Luminosity-Colour 
(PLC) relations (for a more quanti tative discussion of bot h 
random and systematic errors see iBono. fc Cignonil 120051) . 
To illustrate this point, let us assume that the metallicity 
dependence is roughly equal to Mv/S\ogZ ~| 0.4 |. A 
series of random extract ions based on the OGLE cata- 
logue JUdalski et al] 1200 it) indicates that a sample of 1500 
Cepheids with metal abundance ranging from Z=0.002 to 
Z=0.02 and for which i) geometric absolute distances with 
an accuracy better than 2 per cent (~ 0.04 in distance 
modulus); ii) metal abundances with an accuracy better 
than 0.2 dex; Hi) reddening estimates with an accuracy 
better than 0.02 mag are available will allow us not only to 
constrain the metallicity dependence in the optical bands 
with an accuracy better than 0.01 magnitude but also 
to constrain the fine structure of both the PL and PLC 
relations over the entire period range (0.3 < logP < 1.8). 
At the same time, the comparison between theory and 
observations will supply firm constraints on the accuracy 
and plausibility of the physical assumptions adopted in the 
construction of hydrodynamical models of variable stars. 

Classical Cepheids located in the outer disc are impor- 
tant for the distance scale and to estimate disc metallic- 
ity gradients and the shape and scale of Galactic rotation 
JPont et al.lll997l : iMetzger et al.lll99St) . Note that we know 
the chemical composition for only two Cepheid s located at 
distances larger than 16 kpc iLuck et al.ll2003T) and the ra- 
dial velocities for just a handful of Cepheids in the Galactic 
anti-centre. The recent discovery of an extended HI spiral 
arm located in the fourth quadrant of the Milky Way at a 
radial distance of 18 — 24 kpc highlights the potential value 
of a l arge-volume sample of Cephei d tracers in the Galactic 
disc jMcCmre-Grifnths et aljl2004h . 
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Figure 11. Predicted radial velocity and light curves in different 
photometric bands for a Large Magellanic Cloud Bump Cepheid. 
A comparison with empirical I-band data is also shown (third 
panel), taken from OGLE data. 



4-6.2 Low-mass Variable Stars 

RR Lyrae stars are the most popular low-mass distance 
indicators. When compared with Classical Cepheids they 
are at least 2-3 mag fainter but they are ubiquitous in 
the Galactic spheroid and in dwarf galaxies. The optical 
luminosity and metallicity of RR-Lyrae stars are related 
through an My- [Fe/H] relation. Different aspects of this 
relation have been addressed in countless emp i rical and the- 
oretical investigations (see Casteljanjj^j99|; iBond l2003bl : 
ICacciari fc Clementinill2003l : Icatelanl 12001 and references 
therein). The calibration of this relation is still affected by 
uncertainties in distance estimates, reddening corrections 
and evolution off the zero-age horizontal branch (ZAHB). 
However, recent theoretical and empirical evidence indi- 
cates that the K-band Period-Luminosity-Metallicity rela- 
tion is o nly marginally af fe cted by these decepti ve uncer- 
tainties feono et al J 120031) . iBenedict et"aD J2002I) recently 
provided an accurate trigonometric parallax for RR Lyr it- 
self based on observations obtained with the Fine Guide 
Sensor on board HST (although they had only 5 reference 
stars). This notwithstanding, a robust calibration of both 
the zero-point and slope of the My vs [Fe/H] relation still 
awaits better statistics. If we assu me an absolute magni- 
tude of My « 0.6 at [Fe/H] = -1.5 |c acciari fc Clementinil 
I2003T) . we find that the RVS will provide accurate measure- 



ments of chemical composition, reddening, and pulsation 
properties for the entire sample o f metal-poor ([Fe/H] = 
-1.65 ± 0.03: ISuntzeff etafl Il994lh halo RR Lyrae out to 
distances of about 2 — 3 kpc. 

A further interesting aspect for which the RVS will be 
valuable is the detection and measurement of RR Lyrae stars 
in binary systems. We still lack detailed measurements of the 



dynamical masses of RR Lyrae. sinc e at present there are 
only a few Galactic (|Wa de_et a l.ll999T) and Large Magellanic 
Cloud iSoszvnski et al.H2003Tr c"andidates. 

Gaia will also create the unique opportunity to cali- 
brate the PL and the PLC relation for the known Type II 
Cepheids, namely stars of type BL Herculis, W Virginis and 
RV Tauri. At present, the empirical relations are hampered 
by small number statistics because only a han dful of such 
stars have been identifie d in Globular Clust ers dPritzl et alJ 
l2003t) and in the LMC jAlcock et alJll99gT) . These objects 
are very good distance indicators because they are at least 
1 mag brighter than RR Lyrae stars. In addition, theory 
and observations suggest that their propert ies depend only 
marginally on metallicity jBono et alJll997lK 



5 CONCLUSIONS 

The Radial Velocity Spectrometer (RVS) on board the Gaia 
satellite will provide spectra in the 848 nm to 874 nm wave- 
length range with a resolving power R = 11500 for stars 
brighter than magnitude V = 17 over approximately 90 
percent of the sky. The unbiased nature of this data set, 
combined with the multi-epoch observations of each object, 
render the RVS data set unique among spectral surveys of 
the Galaxy. In this paper, we have highlighted the many ar- 
eas in which the RVS spectra constitute an essential comple- 
ment to the astrometric and photometric data which Gaia 
will collect. For studies of Galactic structure and evolution, 
the unbiased sample of radial velocities obtained from the 
RVS spectra will provide the sixth phase-space coordinate 
which is vital for the construction of well-constrained Galac- 
tic models. The multi-epoch nature of the observations will 
allow the identification of large numbers of binary stars in 
the Gaia data set while the availability of spectra taken si- 
multaneously with the photometric observations will be in- 
valuable for the study of detailed stellar properties, including 
stellar variability. 

The Gaia mission, and the RVS in particular, will gener- 
ate a data set of unprecedented size and precision which will 
revolutionise our understanding of all aspects of Galactic as- 
tronomy and stellar astrophysics. Harnessing the full power 
of such a data set presents many challenges, both practi- 
cal and theoretical. For example, as we discussed in § 12.51 
the Hipparcos results for the local velocity distribution em- 
phasise the difficulty of disentangling the global properties 
of the disc distribution function from small scale structure 
such as moving groups. Given this difficulty, it will not nec- 
essarily be fruitful, or even possible, to model the Gaia data 
set via a simple decomposition into known populations such 
as disc, bulge and halo due to the ambiguity in the assign- 
ment of individual stars to specific populations. Instead, a 
global approach will be needed, which makes simultaneous 
use of spatial and kinematic inform ation to mode l the en- 
tire data set as a coherent unit (e.g. lBinnevll2005l) . During 
the coming decade, prior to the satellite launch, much work 
needs to be done in order to determine the iterative Galac- 
tic modelling schemes which will make optimal use of the 
Gaia data set and will lead us to a coherent interpretation 
of the vast quantity of information it will contain. Similar 
developments will be required for the study of binaries and 
stellar astrophysics in the post- Gaia era. 
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